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Abstract 

The splitting processes of bremsstrahlung and pair production in a medium are coherent over 
large distances in the very high energy limit, which leads to a suppression known as the Landau- 
Pomeranchuk-Migdal (LPM) effect. We analyze the case when the coherence lengths of two con¬ 
secutive splitting processes overlap, which is important for understanding corrections to standard 
treatments of the LPM effect in QCD. Previous authors have analyzed this problem in the case of 
overlapping double bremsstrahlung where at least one of the bremsstrahlung gluons is soft. Here 
we show how to generalize to include the case where both splittings are hard. A number of tech¬ 
niques must be developed, and so in this paper we simplify by (i) restricting attention to a subset 
of the interference effects, which we call the “crossed” diagrams, and (ii) working in the large-Ac 
limit. We first develop some general formulas that could in principle be implemented numerically 
(with substantial difficulty). To make more analytic progress, we then focus on the case of a thick, 
homogeneous medium and make the multiple scattering approximation (also known as the q or 
harmonic approximation) appropriate at high energy. We show that the differential rate dT /dx dy 
for overlapping double bremsstrahlung of gluons with momentum fractions x and y can then be 
reduced to the calculation of a 1-dimensional integral, which we perform numerically. [Though this 
paper is unfortunately long, our introduction is enough for getting the gist of the method.] 
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I. INTRODUCTION 


At high enough energy, particles passing through matter (cosmic rays through the at¬ 
mosphere, high-energy partons through a quark-gluon plasma, electrons through an electro¬ 
magnetic calorimeter) lose energy primarily through splitting: hard bremsstrahlung or pair 
production which, when repeated, produces a shower of lower energy particles. Naively, one 
would calculate the rate F for each such splitting by (roughly speaking) computing the cross- 
section a for splitting during a collision between the high energy particle and a particle in 
the medium, as in the bremsstrahlung process depicted in £g. 1, and then finding F ~ nva, 
where n is the density of things to scatter from and v the relative velocity. The flaw in this 
argument, as known since the 1950s, is that the duration of this splitting process, called the 
formation time, grows with energy. At high enough energy, the formation time exceeds the 
mean free time between collisions with the medium, and then consecutive scatterings from 
the medium may no longer be treated as quantum mechanically independent for the purpose 
of computing the splitting rate. This is known as the Landau-Pomeranchuk-Migdal (LPM) 
effect [1, 2], which dramatically suppresses the splitting rate at very high energy. In the 
language of Feynman diagrams, the LPM effect represents important interferences between 
splitting before and after a sequence of elastic collisions with the medium, such as shown in 
£g. 2a for QED and £g. 2b for QCD. The analysis of the LPM effect in QCD was pioneered 
in the 1990s by Baier et ah [3, 4] and Zakharov [5], known collectively as BDMPS-Z. 

A natural question that arises is whether consecutive splittings of the high energy particle, 
and not merely consecutive collisions with the medium, occur within the formation time. 
That is, once we use the LPM or BDMPS-Z formalism to compute the rate for a single 
splitting, can we then treat consecutive splittings as independent and so simply use the 
single-splitting rate in a Monte Carlo to compute the development of the shower (which we 
need for answering detailed questions about energy loss)? Or is there instead a significant 
contribution from processes where the formation times associated with consecutive splittings 
overlap, as depicted in fig. 3? Formally, the probability of a splitting is parametrically of 
order a {as in QCD) in one formation time, and so the formation time is one power of a 
smaller than the typical time between consecutive splittings, as depicted in fig. 4. This means 
that the contribution of overlapping formation times, as in fig. 3, is formally suppressed by 
one factor of a. However, in QCD applications, (i) this Og is at best moderately small at the 
scales relevant for the high-energy splitting process in actual relativistic heavy ion collisions, 
and (ii) this factor of as is accompanied by a potentially large double logarithm [6-8]. It is 
therefore important to analyze the effect of overlapping formation times. 

Note that in the last hgure we have drawn only the high-energy particles. This will be our 
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FIG. 1: High-energy bremsstrahlung (blue) during a collision (black) with a particle from the 
medium. The curly line ending in a cross represents the electromagnetic or gluonic fields in the 
medium created by sources, such as by a nucleus or a passing thermal parton. 
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FIG. 2: Examples of interference terms contributing to the LPM effect in (a) QED and (b) QCD. 
Eor consistency with future hgures, we have used blue for high-energy particles in the amplitude 
and red for high-energy particles in the conjugate amplitude. 


formation times 
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EIG. 3: Two consecutive splittings that are close enough that their formation times overlap. Each 
formation time region is depicted by a green or blue, hatched oval. 


convention throughont the rest of the paper: the collisions with the medium are implicit, and 
the high-energy particle lines should be understood as decorated with numerous interactions 
with the medium, as in £g. 2 . 

Our goal in this paper is to develop formalism for a full computation of the LPM effect 
(to include the treatment of overlapping formation times) in the case of two consecutive 
splittings, such as £g. 3. We will then implement this formalism to compute results for 
double-splitting rates in QCD in certain simplifying limits, to include the multiple-scattering 
approximation (also known as the q or harmonic oscillator approximation), which is appro¬ 
priate for typical events at high energy when the medium is much thicker than the mean 
free path for collisions with the medium. 

Previous authors [ 6 - 8 ] have performed explicit calculations for QCD to leading-log order 
in the limit of y x 1, where x and y are the momentum fractions of two of the hnal 
gluons. [That is, if the initial parton energy is E, the three daughters after the two split¬ 
tings have energies xE, yE, and {l—x—y)E.] These results have interesting and important 
consequences, which we will briefly mention later. Wu [ 8 ] has also developed formalism for 
studying the somewhat more general case of x, ?/ -C 1 without assuming y x, but so far 
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FIG. 4: The hierarchy of scales for typical consecutive splitting, if the relevant a is small and if 
one ignores logarithmic enhancements in QCD associated with one of the daughters becoming soft. 
(The cartoon in this picture assumes that the momenta of the two bremsstrahlung gauge bosons are 
parametrically similar. In QCD, formation times shrink as bremsstrahlung gluons become soft.) 


has only performed explicit calculations when y -C x.^ In this paper, we study the problem 
of general x and y, without requiring that either be small, and we will carry out explicit 
calculations in that case. In the context of the multiple-scattering (harmonic oscillator) 
approximation, we will be able to go beyond leading logarithms and develop methods for 
explicitly computing the full result for double splitting. 


A. What we compute 

In this, our hrst paper on the case of general x and y, we will simplify the discussion by 
focusing on the case of the large Nc limit of QCD. This will simplify treatment of the color 
dynamics of high-energy partons (a simplihcation that is not needed in the y x 1 limit 
considered by previous authors).^ For the sake of a minor, further simplihcation, we will 
also focus on the case where the initial high-energy particle is a gluon rather than a quark 
(and so, in the large-iVc limit, all the daughter particles after each splitting are gluons as 
well). 

Though we will develop formal results that are more general, our explicit calculations in 
this paper will be specialized to the case of a thick, static, homogeneous medium. “Thick” 
means that the medium is wide compared to the typical formation length. 

Finally, in this hrst paper, we will only compute a subset of the contributions to the 
double splitting rate, depicted in hg. 5. It is convenient to also draw these same interference 
terms with the amplitudes and conjugate amplitudes tied together, as in hg. 6. We will 
refer to these contributions as the “crossed” contributions, since they involve two crossed 
lines when drawn as in hg. 6.^ An example of an interference contribution that we will not 


^ There have also been some attempts to include a rough, heuristic treatment of overlapping formation 
times in the Monte Carlo generator JEWEL [9]. We have not yet worked out its relationship to the 
results described here or in refs. [6-8]. 

^ Blaizot and Mehtar-Tani [6] also discuss the simplification afforded by the large-W: limit for moving the 
discussion beyond the soft-gluon limit. 

^ Although these diagrams are crossed when drawn as time-ordered diagrams for the rate, as in fig. 6, 
they are not suppressed in the large-W limit. They are still planar diagrams as far as color factors are 
concerned because one of the crossing lines could have instead been connected by drawing it routed around 
the outside of the oval. (If the solid lines in figs. 5 and 6 were quarks instead of gluons, with x and y 
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+ permutations of (x, 1-x-y ) 

FIG. 5: The subset of interference contributions to double splitting that we will evaluate in this 
paper: the “crossed” diagrams. To simplify the drawing, all particles, including bremsstrahlung 
gluons, are indicated by straight lines. The long-dashed and short-dashed lines are the daughters 
with momentum fractions x and y respectively. The naming of the diagrams indicates the time 
order in which emissions occur in the amplitude and conjugate amplitude. For instance, xyyx 
means first (i) x emission in the amplitude, then (ii) y emission in the conjugate amplitude, then 
(iii) y emission in the amplitude, and then (iv) x emission in the conjugate amplitude. 



+ conjugates 

+ permutations of (x, 1-x-y ) 


FIG. 6: An alternative depiction of fig. 5, with amplitudes (blue) and conjugate amplitudes (red) 
sewn together. The dashed lines are colored according to whether they were first emitted in the 
amplitude or conjugate amplitude. 


calculate in this paper, but plan to address in future work, is the un-crossed diagram of £g. 
7. (The evaluation of this contribution involves some additional subtleties beyond what is 
required for the crossed contributions, to be briefly mentioned in our conclusion.) In this 
paper, we will also, for simplicity, focus on the rate for a high-energy particle to split into 
three on-shell daughters, which all persist after the double-splitting process. That is, we 


being bremsstrahlung gluons, then the diagrams explicitly shown would be 1/Nc suppressed, but some of 
the “permutations” referred to in the figure, such as the one shown later in fig. 9, would not be.) 
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FIG. 7: An example of an interference diagram not evaluated in this paper, to be treated in later 
work. 



FIG. 8: An example of a virtual loop correction to single splitting. 


will not yet consider loop corrections to single splitting, such as shown in hg. 8 . Such loop 
corrections are necessary for using calculations of splitting rates to compute energy loss, as 
has been analyzed in refs. [ 6 - 8 ] for the case x, ?/ -C 1 . 

The “+ permutations” line of hgs. 5 and 6 contains more types of interference terms than 
one might at hrst realize. For example, it contains the diagram of £g. 9 corresponding to 
the interference of (i) separate emission of x and y with (ii) the emission of an x+y gluon 
that then splits into x and y. This diagram is equivalent to the x -H- 1—x—y permutation 
of the xyyx interference shown explicitly in hg. 5. We will correspondingly refer to hg. 9 as 
the zyyz interference, where z = \—x—y. 


B. Overview of method and what’s new 

Here, in broad outline, are the elements involved in the calculation. We start by turning 
the problem into a problem in two-dimensional non-Hermitian quantum mechanics, slightly 
generalizing techniques that have been used by other authors. 


1. 3- and 4-po-Tzticle quantum mechanics 

Start by considering an interference contribution to the rate for single splitting, shown 
in hg. 10. In the time interval < t < between splitting in the amplitude and conjugate 
amplitude, the number of high-energy particles in the diagram does not change. As we shall 
review later, at high energy it is possible to reduce the problem in this time interval to time 
evolution in two-dimensional quantum mechanics for three non-relativistic particles, two 
representing the two daughter particles in the amplitude in hg. 10 during this time interval 
and one representing the parent particle in the conjugate amplitude in hg. 10 during the 
same time interval. The two dimensions of the quantum mechanics problem are the plane 
transverse to an axis approximately aligned with the initial direction of the parent. The 
“non-relativistic” character of the transverse dynamics comes from the large-p^ (and small- 
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FIG. 9: One way to draw the zyyz interference contribution, which is the x 
permutation of the xyyx interference shown explicitly in fig. 5. 


= l-x-y 



FIG. 10: The basic interference contribution for 
labeling convention for the splitting times in the 



single (independent) splitting, together with our 
amplitude and conjugate amplitude. 


angle) approximation to the kinetic energy Sp of a high-energy particle of mass Ai, 

= \ pI+P± + —Pz + -— 77 ^ + constant. (1.1) 

V 2p^ 

The right-hand side looks like a two-dimensional non-relativistic kinetic energy p\/2m with 
a “mass” m with magnitude |pz|, and (for our purposes) \pz\ may be treated as constant 
in the above approximation.'^ As we’ll discuss, the 3-particle quantum mechanics problem 
that describes fig. 10 has the funny feature that the “mass” m which describes the particle 
in the conjugate amplitude is taken to be negative, which arises from the sign difference 
between the evolution exp(—iift) of amplitudes vs. the evolution eyip{+iHt) of conjugate 
amplitudes. Because of overall longitudinal momentum conservation in the splitting process, 
the two positive and one negative “mass” in this quantum-mechanics problem satisfy the 
unusual property (from the perspective of non-relativistic quantum mechanics) that 

mi + m 2 -\r m 3 = 0 . ( 1 - 2 ) 

Another unusual feature of the effective quantum mechanics problem is that it has a non- 
Hermitian Hamiltonian. As the particles propagate, they interact with the background helds 
in the plasma. For a thermal medium (or other random medium), these background helds 
are random with well-dehned correlations. When computing a rate in a random background, 
one should average the rate over that randomness (e.g. average over the thermal ensemble). 
After this average, the interactions of the particles with the background helds, and their 


^ The distinction “for our purposes” is important. In gp ~ pz +p\/2pz the variation of the first term 

Pz can be as large as the second term p\/2pz in our context. However, the leading pz term will cancel in 
the combinations of Cp’s that will appear in our problem, such as (2.14). 
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FIG. 11: The xyyx interference contribution of figs. 


5 and 6 with individual splitting times labeled. 


correlations, may be replaced by a “potential energy” term in the qnantum mechanics prob¬ 
lem, as we will later review more concretely.^ However, unlike a normal quantum mechanics 
problem, this potential energy is not real-valued. The effective quantum mechanics problem 
that reproduces the medium-averaged splitting rate (rather than the unaveraged splitting 
amplitude) has a non-Hermitian Hamiltonian. The non-Hermiticity of the Hamiltonian ac¬ 
counts for quantum decoherence over long times: it causes the interference contribution of 
£g. 10 to decay exponentially for time separations t^—t^ large compared to the formation 
time. 

Now consider a particular interference contribution to the double splitting rate, such as 
the xyxy interference shown in fig. 11. Within each time interval between splittings in the 
figure, such as (a) t^ < t < ty, (b) ty < t < ty, or (c) ty < t < t^, the number of high- 
energy partons does not change. Using the same methods outlined above, we can reduce the 
problem of medium-averaged time evolution during these intervals to two-dimensional non- 
Hermitian non-relativistic quantum mechanics for (a) three, (b) four, and (c) three particles 
respectively. For the 4-particle evolution for ty < t < ty, the “masses” of the quantum 
mechanics problem will satisfy 


mi + 1712 + 'm 3 + ^4 = 0, 


(1.3) 


in analogy to the 3-particle case (1.2). 

If we can solve the time evolution in the non-Hermitian quantum mechanics problems, 
we can glue the results together using the field-theory matrix elements for the splittings at 
t = tx, ty, ty, and tx in £g. 11 to compute the result that we want for that interference 
contribution. We may similarly compute all the other interference contributions and add 
them together. Fortunately, for the sake of making the calculation practical, there is a major 
simplihcation to the problem, which we now outline. 


® A potential in non-relativistic quantum mechanics represents an instantaneous correlation in our problem. 
This is a valid approximation in our problem because the time it takes for high-energy particles to cross 
one correlation length of the medium is parametrically small compared to the formation time in the high- 
energy limit. The dynamics of the high-energy particles on the time scale of the formation time is what 
is relevant to our calculation. 
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2. Simplification to 1- and 2-particle quantum mechanics 


We will need to solve for two-dimensional quantum mechanical evolution in problems of 
the form 

where h* are the two-dimensional transverse positions of the particles, Pi are the correspond¬ 
ing transverse momenta, and V is the (non-real) potential that incorporates the medium- 
averaged effects of scattering from the medium, which may depend on t ^ z (where 2 : in 
this context refers to the longitudinal coordinate 0 : 3 ). At high energies, the medium will 
be (statistically) translation invariant in the transverse plane over the small transverse dis¬ 
tances relevant to the problem, and so the potential V will be similarly invariant. For generic 
masses nii, that invariance lets us factor out the analog of the trivial “center of mass” motion 
of the problem (which turns out not to be relevant) and so reduce the Wparticle problem 
to an effective A^—1 particle problem. However, as we shall discuss, something very special 
happens when ~ 0 ) which will allow us to eliminate yet another two-dimensional 

degree of freedom and so reduce the N particle problem to an N—2 particle problem. So, 
3-particle evolution can be reduced to a simple problem of 1-particle quantum mechanics, 
which is equivalent to simplihcations originally used by BDMPS-Z to analyze the problem 
of the single-splitting rate via £g. 10 in the case of general x. In our case, it will allow us (at 
the expense of having to develop some formalism) to reduce the 4-particle two-dimensional 
quantum mechanics problem required for the ty < t < ty region of fig. 11 to a simpler 
2 -particle two-dimensional quantum mechanics problem. 


3. Further specialization 

In the multiple scattering (harmonic) approximation, appropriate at high energy both 
for thick media and for “typical” events in relatively thin media,® the potential V in (1.4) 
becomes a (non-Hermitian) harmonic oscillator potential. Then we will just need to (i) 
solve for time-evolution of 1 -particle and 2 -particle harmonic oscillator problems; (ii) glue 
the results together with splitting matrix elements; (hi) integrate over the four times ty, 
tx, and ty associated with those splittings in £g. 11 , for example; and (iv) integrate over 
all relevant intermediate states (i.e. transverse positions) of the high-energy particles at the 
intermediate times ty and ty in fig. 11. That sounds like a lot of complicated integration. 
However, for the case of a thick, homogeneous medium, we will see that it is possible to do 
all of these integrations in closed form, except for a single time integral over the duration 


® Here, thin and thick are defined relative to the typical formation length in an infinite medium, which 
for QCD is ~ i/w/g, where to is the smallest energy of the daughters. There is often a great deal of 
confusion about the applicability of the harmonic oscillator approximation to thin media because interest 
is often focused on calculations of the medium effect on the average energy loss (AE ), which in the thin- 
media case can be dominated by rare events with relatively large AE. However, the harmonic oscillator 
approximation still describes typical events in the thin-media case, provided the thickness of the medium 
is large compared to the mean free path for collisions. See the notes and references given for (2.22) in 
appendix A. 
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At of the 4-particle evolution (i.e. over At = ty — ty for fig. 11). The hnal At integral will 
be performed numerically. 


4- Other aspects worth noting before we begin 


Before we launch into the main body of the paper, we will hrst note some other aspects 
of our problem (double splitting for general x and y) that did not arise in previous work 
(specihc to small x and y). The first is the necessity for a careful treatment of the helicity 
of high-energy particles between splittings. In the calculation of the single-splitting rate of 
fig. 10, it is possible to express the x dependence of the splitting matrix elements at and 
tx in terms of spin-averaged Dokshitzer-Gribov-Lipatov-Altarelli-Parisi (DGLAP) splitting 
functions. As an example, the differential splitting rate in the case of a thick, homogeneous 
medium is [10] 

^ P{x) \no\, (1.5a) 

dx Try 2 

where P{x) is the relevant spin-averaged DGLAP splitting function and Dq is the (complex) 
frequency associated with the non-Hermitian harmonic oscillator in the multiple scattering 
approximation. For g ^ gg splitting for example, Do turns out to be 


Do — 



(1.5b) 


{qr, known as the jet-quenching parameter, is the average squared transverse momentum 
per unit length that is transferred from the medinm to a high-energy particle with color 
representation R.) When we stndy double splitting, snch as in hgs. 5 and 6, the helicities 
associated with the various splittings are tied together in a non-trivial way that cannot be 
packaged into spin-averaged DGLAP splitting fnnctions for general valnes of x and y. One 
of our tasks will be to distinguish different intermediate-state helicity amplitudes in our 
treatment of the splitting matrix elements. 

The other non-trivial feature that we want to mention is that individnal interference con- 
tribntions, such as each of those in £g. 5, will tnrn ont to have an nnphysical logarithmic 
UV divergence arising from coincidence of three of the four splitting times in the diagram.^ 
Since the diagrams for the amplitnde are tree level (and similarly for the conjngate ampli- 
tnde), there are no loops involved in the amplitnde, and so there shonld be no actnal UV 
divergences in the final result for the differential donble splitting rate dT/dx dy once we snm 
over all contribntions. We will see that the short-time UV issne is resolved by being carefnl 
abont operator ordering and ie prescriptions. 


C. What we find 

Since we are only compnting the “crossed” snbset of contribntions to the donble-splitting 
rate (fig. 6) and not yet all contribntions, the results of this paper lie in the development of 


This is not the / d{At)/At divergence that is part of the double log behavior found for y <C x in refs. 
[6-8]. For fixed x and y (which is what we consider in this paper), that particular logarithm is cut off at 
small At ~ y^yE/xq, as we review in section IX A. 
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calculational methods rather than a hnal, total answer for the differential double-splitting 
rate dV/dxdy. That said, it is still interesting to ask if the contribution of this subset to 
the result for dT/dxdy has any interesting qualitative features. First, we will check that for 
1 / a; 1, we hnd a contribution with y dependence of the form 


dr 

dx dy 


oc-h 


J crossed 


Qs Mv/x) 

X^l‘^y 


+ ••• . 


( 1 . 6 ) 


When integrated over y, this gives a double log correction oc o;^ log^. We will see (with 
some caveats explained later) that this result matches recent y ^ x results by Blaizot and 
Mehtar-Tani [6] and lancu [7] and Wu [8], after accounting for the fact that we have not 
included virtual corrections like £g. 8. The x~^f‘^y~^ \\i{%)/x) contribution to (1.6) arises from 
interference diagrams of the form xyyx, xyyx, zyyz and zyyz, where the zyyz interference 
was depicted in fig. 9. As discussed in refs. [6-8], this correction is due to small-a; physics in 
the collisions of the high-energy particles with the medium, and it may be absorbed into a 
running of the parameter q that characterizes the strength of interactions of the high-energy 
particles with the medium in the multiple-scattering approximation [6, 11, 12]. 

However, in the limit of small y and fixed x, we will find an even more divergent contri¬ 
bution from crossed diagrams that scales with y as 


dT ^ Q^s Mv/x) 
dxdy\^^^^^^^ xy^/^ 


(1.7) 


The physical meaning of this soft divergence will be discussed in future work. 


D. Outline 

At the end of this introduction, we will briefly comment on the relevant scale of the ag that 
controls a perturbative treatment in the number of overlapping hard splittings. In section II, 
we review formalism and results for computing single splitting processes, as in fig. 10, which 
will set up some of the language and formalism we later need for double splitting. Section 
III turns to how to use symmetry to reduce the problem of N particle quantum evolution 
in interference diagrams to an effective quantum mechanics problem involving only N—2 
particles. We apply that technique to the double splitting problem in section IV, which also 
discusses the rules of diagrammatic perturbation theory in this language, how the large-A^,. 
limit simplifies the discussion of color in this paper, and the relationship of splitting vertex 
matrix elements to helicity-dependent DGLAP structure functions. In that section, we go 
as far as we can without either attempting complicated numerics or introducing additional 
approximations. To go further, we turn to the multiple scattering (harmonic oscillator) and 
thick-media approximations in section V, where, as a first example, we reduce the calculation 
of the xyyx interference in fig. 11 to the form of a single integral over At = ty — ty. We 
also show that this integral is divergent as At —)■ 0 and isolate the divergence. In section 
VI, we show how all the other “crossed” interference diagrams of fig. 6 can be related to 
the xyyx result, and we show that the small-At divergences cancel among these diagrams. 
We then discuss in section VII that the canceling divergences nevertheless leave behind a 
hnite piece associated with the residue of a pole at At = 0. Getting this right requires a 
discussion of the appropriate ie prescriptions associated with At —)■ 0. Section VIII gives 
a summary of our final results in terms of a single convergent At integral which can be 
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performed numerically. In section IX, we discuss our results in the limit where at least one 
of the radiated gluons becomes soft and show agreement between the pieces of our results 
that match up with pieces of previous authors’ soft-gluon results. Finally, section X offers 
our slightly redundant but extremely short conclusion. A number of details are left for 
appendices. 


E. Note on the scale of Qg 

The size of ttg depends on momentum scale, and there are two parametrically different 
scales relevant to the size of couplings in splitting processes such as £g. 2. One is the 
strength of the interaction of high-energy particles with the medium (the couplings to the 
black, vertical gluon lines in fig. 2). The smallest momentum transfers that are important are 
of order the inverse Debye screening length, mo, which is a characteristic of the medium. In 
theoretical studies, one may or may not want to consider as{mY)) as small, but it is certainly 
not very small in quark-gluon plasmas created in experiment. The other important scale is 
the one associated with the high-energy splitting vertex in fig. 2. The relevant scale there is 
the typical transverse momentum Q± between the daughters, which depends on the energy. 
(This is the scale of the momenta of each daughter in the daughters’ center-of-momentum 
frame.) In the simple case of hard bremsstrahlung in a thick medium, the LPM effect 
causes Q± to grow slowly with energy as® Q± ~ {qEY^*. If E is large enough, the q;s(Q_l) 
characterizing the splitting vertex may therefore be small, regardless of whether aY^D) 
is small. This approximation—that the amplitude for each high-energy splitting vertex is 
small—is one that we will use implicitly throughout this paper. The factors of as in this 
paper will all implicitly be factors of q;s(Q±) unless explicitly noted otherwise. 


II. REVIEW OF SINGLE SPLITTING 

We will set the stage by summarizing, in our own language, the starting point for the 
basic formalism (used in one form or another since BDMPS-Z) describing the LPM effect 
in QCD for single splitting as in fig. 2. Everything in this section is equivalent to results 
already known in the combined literature on this problem, but it will be useful for hxing 
our own conventions and in particular fixing the language we will use when we move on to 
discuss double splitting.® 


® The relative Q± of the transverse momentum between the daughters is of order the total momentum 
transfer from the medium during the formation time tform: and so Q\ ~ (?horm- For hard (x ~ 1 ~ 1 — cc) 
bremsstrahlung in a thick medium, the LPM effect gives tform ~ i/w/g ~ i/E/g, and so Q± ~ {qEY^'^. 

® Our formulation of the problem owes the most to Zakharov [5], together with elements from refs. [2, 13, 14]. 
The notation we use is closest to ref. [14], whose appendix includes a translation table to the notation of 
several other works. Other noteworthy formalisms for studying LPM in QCD include BDMPS [3, 4], but 
also, in various limits or with various extensions, GLV [15], ASW [16], and higher twist [17]. See ref. [18] 
for comparative discussion and further references. 
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FIG. 12: Our conventions for labeling momenta for the interference contribution for single splitting. 
In the QED case, fcx = fcf since the photon does not directly interact with the medium. 


A. Relating bremsstrahlung to 3-particle non-Hermitian quantum mechanics 


Let’s first discuss the case of QED, as in fig. 12 with the dashed line representing a 
bremsstrahlung photon. Let 5H be the part of the Hamiltonian that contains the splitting 
vertices for the high-energy particles. Then, working to leading order in 5H, the differential 
bremsstrahlung probability dl/dx for producing a photon with momentum fraction x is^° 


dl 

dx 



(Pffcf, 4|Pxfcf, 4) (Pxfcf I dH\pi) 


'Px 


X 


(Pf fcf I -z (5iL IPx) (px, 4 |pi, 


'Px 



( 2 . 1 ) 


where the p is the transverse momentum of the charged particle and k is the transverse 
momentum of the bremsstrahlung photon, labeled as in Eg. 12. The high-energy particles 
propagate through gauge helds sourced by the medium, and the double angle brackets denote 
statistical averaging over those gauge helds. The factors in the hrst parenthesis give the 
amplitude in £g. 12, which consists of a splitting matrix element (px^f|—z|Pi) at the 
hrst time G, followed by the propagation (pffcf, txipxfef, G) in the background held from G 
to tx. In QED, the bremsstrahlung photon (fc) is just a spectator during this propagation. 
The factors in the second parenthesis give the conjugate amplitude in hg. 12, which has 
propagation hrst, followed by splitting at the later time t^- The overall 2Re{- • • } combines 
(i) the interference term of hg. 12, which is the case where emission happens hrst in the 


In this paper, we do not consider bremsstrahlung associated with whatever earlier hard process created our 
initial high energy parton in the medium in the first place. We consider only late enough times (t ^ 1 /rei 
where Tei is the mean free path for collisions with the medium) that the high-energy partons are already 
approximately on-shell. (By the uncertainty principle, a particle having collisions roughly every time 
interval Tei will be off-shell in energy by AE ^ l/rei, which corresponds to p^Pfi E AE ^ Ejr^x^ 
which is what we mean here by “approximately on-shell.”) For a recent interesting discussion of the 
interplay between early vacuum bremsstrahlung and late medium-induced bremsstrahlung, see Kurkela 
and Wiedemann [19]. 
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amplitude, with (ii) its complex conjugate, which is the case where emission happens hrst in 
the conjugate amplitude. The sum in (2.1) is over hnal state polarizations of the daughters. 
The symbol V± represents the (inhnite) 2-volume (area) of the transverse plane, upon which 
hnal results will not depend. 

In this paper, we will delegate additional discussion of some formulas to appendix A. In 
the case at hand, see appendix A if wondering why one need not include time evolution of 
the hnal state \pfkf) after in (2.1), nor time evolution of the initial state |pi) before tx- 
We also comment there on the prefactor E/27iV± in (2.1). 

It is useful to rewrite the conjugate amplitude to express (2.1) as 


dl 

dx 



(Pffcf, txIPx^f, fx) {p^h I -i dH\p-^ 


'Px 


X 


(Pi, 4|Px, 4)(PxK dH\piki) 


'Px 


( 2 . 2 ) 


Time-evolution in the amplitude is given by the factor 

{pffcf,fx|PxA;f,tx) = (^2.3) 

where H( 2 ) is the Hamiltonian for propagating two high-energy particles through the 
medium; the charged particle and the photon. (Terms 5H in the full Hamiltonian which 
change the number of high-energy particles are not part of H{ 2 )-) The time evolution in the 
conjugate amplitude is analogously in the factor 

(Pi,tx|Px,tx) = (2.4) 


where if(i) is the Hamiltonian for propagating one high-energy particle through the medium. 
Since the interactions with the medium appear only in and H{ 2 ), the medinm-average 
in (2.2) can be restricted to the pieces 


— (ix—ix) 


Pxki){Pi\ 




(2.5) 


This is just time evolution of an initial \Pxkf){p[\ starting from time t^- The object \Pxkf){Pi\ 
lives in the space 

He «) He,7 = He (g) He (g) H.^, (2.6) 

where Hg is the Hilbert space of states of a high-energy electron, and He,.y is the Hilbert space 
of states with both a high-energy electron and high-energy photon. Instead of thinking of 
(2.6) as a product of Hilbert spaces, it is convenient to formally think of a single Hilbert space 
of three particles: one electron, one photon, and one conjugated electron.Correspondingly, 
we rewrite (2.5) in the form of a 3-particle evolution 

\pi,P2, P3) )), (2.7) 

where 77(1+2) = 77(2) — -^(i)- Here 77(2) operates on the two particles associated with the 
amplitude and 77(i) on the one associated with the conjugate amplitude. It will be convenient 


This picture is at the heart of Zakharov’s early work [5]. 
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to choose the convention that a Pi above that represents the momentum of a particle in the 
conjugate amplitude has been negated. So, for instance, (2.5) may be recast in the form 
(2.7) with the identideation 

(Pi,P2,P3) = (-Pi,Px, k) aXt = 4- (2.8) 

(Our convention is that the “first” of the three particles is the one from the conjugate 
amplitude.) With this sign convention for the Pi, momentum conservation Pi = Px + kf in 
fig. 12 means that 

Pi+P2+P3 = 0. (2.9) 

The medium average in (2.7) only affects the evolution operator (which depends on the 
background fields), and so we may rewrite (2.7) as 

e-*^(*""*")|pi,P2,P3) = ((e"*^(J+2)^‘""*"^^^|pi,P2,P3), (2.10) 

where "H is going to be the Hamiltonian of our effective 3-particle quantum mechanics 
problem. Note from this definition that "H need not be Hermitian, even though i7(i+2) is. 


B. Brief review of the form of H 


As we’ll very briefly review, the form of "H turns out to be (1.4), 

with non-relativistic “masses” given in terms of the initial parton energy E by 

rui = XiE. (2.11b) 

We have used a similar minus sign convention for defining longitudinal momentum fractions 
Xi as we did above for defining transverse momenta pp momentum fractions Xi appearing 
in the conjugate amplitude are negated. In our case here, 

{xi,X 2 , X 3 ) = (-1, 1 -x, x), (2.12) 


and so mi is negative while m 2 and m 3 are positive. Note that, with this sign convention, 

Xi + X2 + X3 = 0, (2.13) 

which is equivalent to the mass relation (1.2). 

The kinetic terms in (2.11) can be understood by taking the limiting case where there is 
no medium, so that the high-energy particles are free. Then 




( 2 ) 


kf(i) — (£p2 + ^pa) — £; 


Pi 5 


( 2 . 14 ) 


where £p is the energy of a single free particle. Using the large-p^ expansion (1.1) in (2.14) 
yields the kinetic term of (2.11a) plus a term 


Mf 

2xiE ' 2X2E ' 2X3E 


+ 


, Ml 


+ 


( 2 . 15 ) 
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for the potential. In this paper, we will assume for simplicity that the energy is large 
enough that effects of parton masses Mi (including effective parton masses in the medium) 
are negligible. 

The dominant effect of interactions with the medium is to add an imaginary part —iT^ 
to l-i related to the differential rate of scattering from the medium. For the case of a 
weakly-coupled, QED plasma, this turns out to be 

V (hi, hs, hs; t) = -i [e^ rei(0, t) - rei(h2-hi, t)] , (2.16) 

where e is the charge of the charged high-energy particle, corresponds to the photon. 
Above, Fei = e^rei(O) is the rate of elastic scattering from the medium, which will vary 
with time if the medium properties along the path of the high-energy particles vary with 
time. The generalization rei(h, t) is dehned as the Fourier transform of the differential rate 
of scattering dVei/d^qi. with respect to the transverse momentum transfer q±. 

r.,(6) . / ,v ^ i/dV (2.17) 

The second term in (2.16) corresponds to background held correlations between amplitude 
and conjugate amplitude in £g. 2a. The first term corresponds to self-energies of the charged 
particle lines arising from correlations between the amplitude and itself, or between the 
conjugate amplitude and itself.The relative sign in (2.16) arises because the second term 
corresponds, in the language of "H, to the interaction of a charged particle (charge e) and a 
conjugated charged particle (charge —e). 

To generalize to strongly-coupled plasmas (but weakly coupled splitting matrix elements 
5H for high-energy particles), it has been argued that the potential V may instead be 
related to the value of real-time Wilson loops which contain two long, parallel, light-like 
lines separated by 6 = &2 — hi [21]. In this paper, we will not commit to weakly or strongly 
coupled plasmas. We will keep V general in the first part of our analysis, before later 
specializing to the multiple scattering approximation. 


C. Generalization to QCD 

The case of QCD is similar, but the bremsstrahlung gluon also interacts with the medium, 
and so correlations also arise between the bremsstrahlung gluon (bs) and the other two 
particles (bi,b 2 ). For a weakly-coupled plasma, the corresponding potential in (2.11a) is 


For a relativistic plasma and for zero-temperature masses small compared to T, this condition is formally 
that the energies be 3> T. But mass effects can be important for thin media in the harmonic oscillator 
approximation [20]. 

The correlation of two interactions in the amplitude has a real part in addition to the imaginary part, 
but this may be absorbed by replacing the mass of the high-energy particle by its effective mass in the 
medium. 
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V{b,, b2, bs; t) = -t + T| + T|) fei(0, t) + ■ T, f,,{b2-b,,t) 

+ • T 2 fel(b3 —^ 2 , t) + (yf^Ti ■ T 3 fel(bl—bs, t) , (2.18) 

where T“ are the color generators that operate on the color space of particle i, and Ti ■ Tj = 
T“Tj^. The differential rate of elastic scattering from the medium for a particle with color 
representation R is normalized here (in weak coupling) as T/j^ei = g'^CnTei. The hrst term 
in (2.18) is the sum of the self-energies of the three particles. The R ■ Tj terms involve 
correlations between particle i and j, with charges gR and gTj respectively. This means that 
the T 2 ■ Ti and T 3 ■ Ti terms correspond to correlations between interactions in the amplitude 
and conjugate amplitude, according to our convention that particle 1 is the particle in the 
conjugate amplitude. The T 2 ■ T 3 involves correlations between background interactions of 
the two daughter particles in the amplitude.^® Color conservation implies that 

Ti+T 2 + T 3 = 0, (2.19) 

similar to pi -|- p 2 + Ps = 0 and Xi + X 2 + x^, = 0. Relation (2.19) can be used^^ to express 
Ti ■ Tj in terms of hxed Casimirs C* = T^, putting (2.18) into the form 

R(bi,b2,b3;f) = -tg^ i(C'i+C' 2 -C' 3 ) Afel(b 2 -bi,t) 

+ ^(^* 2 +^ 3 —Cl) Arei(b 3 —b 2 , t)-|-|(C' 3 -|-C'i— 6 * 2 ) Arei(bi—bs,f) , ( 2 . 20 a) 

where 

Afel(b) = rel(O) - fel(b) = J (1 - ( 2 . 20 b) 

Note that the potential V depends only on differences bi—bj, which is a consequence of 
transverse translation invariance. 

It is worth mentioning in passing a simplihcation in the case of small x: the masses mi 
and m 2 in ( 2 . 11 ) then have large magnitudes compared to m 3 . As a result, there is little 
(transverse) motion of particles 1 and 2 , and so bi and b 2 stay close together compared 
to b 3 . Then particle 2 and conjugate particle 1 can be thought of as forming a tiny, hxed 
color dipole which interacts with particle 3 (e.g. a qq or gg dipole interacting with a softer 
gluon). This is one way that the 3-particle problem can be reduced to a 1-particle problem 
(of the small-x particle). However, we will be considering generic x and will not make such 
approximations in this work. 

The here in section IIC characterizes interactions with the medium, and so the relevant scale for this 
running coupling can be as low as the inverse Debye screening length. (We promised earlier to alert 
readers to those few cases where the strong coupling in this paper was not evaluated at a scale Q± that 
grows with energy E.) 

Note that (gTi,gr 2 , ( 7 T 3 ) —>■ (—e,-be, 0) turns the QCD result (2.18) into the QED result (2.16). Here, 
gTi —>• 0 represents the charge of the photon. 

The effects of the imaginary part of this correlation dominate over those of the real part. 

For example, (2.19) implies ('ri-l-T 2 )^ = T|, which gives ri-T 2 = — |(T^-bT| — r|) = — ^(Ci-b (72 — Ca). 
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D. The (optional) harmonic oscillator approximation 

There is no reason that one cannot work hard to numerically solve the quantum mechanics 
problem with the full potential V [22, 23]. However, at high energy (for a medium thick 
compared to the mean free path for elastic scattering) there is an oft-used simplihcation; it is 
difficult to deflect high-energy particles, and so (given that the particles start out together), 
the values of hi—bj will be small during the formation time. One might therefore expect 
that one can replace the potential V above by a quadratic approximation around bi—bj = 0. 
Formally, the small-b limit of (2.20) gives 


V{bi,b2,b3;t) 


(gi + q2- q3){b2-bif 


+ (g 2 + ga - qi) {b3-b2f + (gs + gi - g 2 ) (hi-hs)^ , 


( 2 . 21 ) 


where qn, the average squared transverse momentum transferred from the medium per unit 
length, is 

quit) = j d?q± ( 2 - 22 ) 

In the case of a strongly-interacting plasma, the form of V may be more complicated 
than the weakly-interacting form of (2.20), but we give an argument in Appendix A that 
(2.21) holds for any quadratic approximation to the small bi—bj limit. 

See appendix A for caveats and further discussion of the dehnition (2.22) of g/j. 


E. Reduction from 3-particle to 1-particle quantum mechanics 

In our application, Pi + P 2 + P 3 = 0, and so there are seemingly only two independent 
momenta in the Hamiltonian "H (2.11a) in the case of interest. However, we get a further 
simplihcation from the 3-dimensional rotation invariance of our problem. Specihcally, we 
could choose our axis to be in a slightly different direction, while still retaining the large-p^ 
approximations that we have made. In the context of those approximations, the answer for 
the splitting rate should accordingly be invariant under transformations of the form 


{Pi.Viz) iPi+Vizi, Viz), (2.23) 

where ^ is an arbitrary vector in the transverse plane. The transformation (2.23) corresponds 
to a 3-dimensional rotation by 0 = x ^ once one linearizes the rotation in \ 6 \. (Here, 
is the unit vector in the direction.) Writing pi^ as rrij = XiE, the above transformation is 

Pi^Pi + rriii =Pi + XiEi. (2.24) 

(In the language of two-dimensional quantum mechanics, this is Galilean boost invariance, 
with ^ playing the role of the relative velocity of the frames.^®) Because the rate we are 
computing is independent of the rotations implemented by (2.24), one expects that the 


Three-dimensional rotations mix (b,z), whereas two-dimensional Galilean boost invariance mixes 
The two correspond to each other here because z ~ t in the large-pz approximation. 
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calculation only depends on invariant combinations of the Pi, such as XiPj — xjPi. (We will 
give a more precise argument later.) Because Pi + P 2 + Ps = 0 and xi + X 2 + x^ = Q m. our 
problem, there is only one independent such combination: 


P = X 2 P 1 - X 1 P 2 = X 3 P 2 - X 2 P 3 = XiPs - X 3 P 1 . 


(2.25) 


And, indeed, the 3-particle kinetic energy term in "H (2.11a) can be rewritten in terms of P 
in the form of a 1 -particle kinetic energy: 

^ + ^ + A = 4 (2.26) 

2mi 2m2 2m3 2M 

with 

M = —X 1 X 2 X 3 E. (2.27) 

From (2.12) M = x{l — x)E, which is positive since our xi is negative. 

The most natural choice for a variable conjugate to XjPi — XiPj is (hj — bj)/{xi + Xj), 
noting that the commutator 


XjPi - XiPj 


bj - bj 

(Xi +Xj)- 


(2.28) 


Given the string of equalities in (2.25), one might expect that in our problem there is 
correspondingly a single independent transverse position variable B with 


bi - b2 _ b2 - bs _ bs - bi 
{X 1 +X 2 ) {X 2 + X 3 ) (a;3-l-a;i)’ 


(2.29) 


This is indeed the case, but we will hold off on giving a more precise argument until section 
III, where we will simultaneously work out the generalizations that we need for the double 
splitting calculation (where we will additionally reduce 4-particle quantum mechanics to 
2-particle quantum mechanics). Note that we could use xi -|- 0:2 - 1 - 0:3 = 0 to simplify xi + X 2 
to —X 3 and so on in (2.29), but the combinations in (2.29) turn out to be the ones that 
generalize nicely to the case of more than 3 particles. 

Let’s now slightly change notation for the translation-invariant 3-particle potential in "H 
by writing 

V (hi, b2, bs; t) = U (hi-h2, h2-h3, b^-bi] t). (2.30) 

Listing three bi—bj combinations on the right-hand side is redundant, since only two are 
independent, but the redundancy makes manifest the permutation symmetries. Taking 
(2.29) at face value for the moment, the 3-particle Hamiltonian (2.11a) then reduces to the 
1-particle Hamiltonian 

p2 

■H = — + U{{Xi+X 2 )B, {x 2 +X 3 )B, {x 3 +Xi)B-,t). (2.31) 


As an example, the potential (2.20a) for weakly-coupled plasmas would give (now using 
xi+ X 2 = -X 3 , etc.)^^ 


n = 


2M 


-w 


|(C'i+C'2-C'3)Arei(-a;3B,t) 

+ \{C2+C3-Ci) Afei(-a;iB,f) + l(C'3+C'i-C'2) Afei(-a;2B, t) 


, (2.32) 


Here again, as in footnote 14, the characterizes interactions with the medium. 
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which is equivalent (with different notation and slight generalization) to the formalism orig¬ 
inally used by BDMPS-Z. As another example (not restricted to weak coupling), the general 
harmonic oscillator approximation of (2.21) gives 


with 



+ B^ 

(2.33a) 


X 2 X 3 J ' 

(2.33b) 


F. Finishing np 


In QCD, bremsstrahlung gluons can interact directly with the medium, and so the inter¬ 
mediate and final momenta and fcf need not be the same in fig. 12. The single splitting 
probability (2.2) for QED correspondingly generalizes to 



(The discussion in this section applies equally well to q ^ qg, g ^ gg, and g qq, though 
in this paper we will eventually focus solely on gluons.) We will later explicitly discuss the 
details of reformulating this formula in terms of P while getting all the factors right. For 
now, we jump ahead to review the general expression for the final result, which is 


(£_ 

dx 


aPi^sjx) 

[x{l-x)E]‘^ 



P (P- t-\P t ) 

-*• X X5 ^X|-*- X? ^x/5 


(2.35) 


where Pi^ 3 {x) is the usual (spin-averaged) DGLAP splitting function. In this formula, 
(Px, tx|Px, tx) is the propagator of the 1-particle quantum mechanics problem of (2.31). The 
two SH matrix elements in (2.34) are responsible for the two factors of P in (2.35) and for 
the factor of Q;Pi^. 3 (a;). One may alternatively work in B space instead of P space, in which 
case the above formula is [5] 


dl aPi^s^x) 


dx [x{l — x)E]‘ 


Re 


dp dtx ^Bx ' G) 




Bx=Bx=0 


(2.36) 


We will find it convenient to work in B space when later deriving our results for double 
splitting. In the case of a thick and homogeneous medium, (2.36) reduces to 


dT aPi^ 3 {x) 


dx [a;(l — x)E]‘^ 


roo 

Re / d{At) Vbx ■ Vb^(Px, At\B^, 0) 
do 


Bx=Bx=0 


(2.37) 


for the differential rate of single splitting. 
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In order to make use of this formula, one could calculate numerical results for the prop¬ 
agator [22, 23]. With the harmonic oscillator approximation, however, one can be more 
analytic. For example, for a thick, homogeneous medium, 

{B, At\B', 0) = - ° ^ ° ^ exp^|Mffo [{B"^ + B'"^) cot(f2o^^) — 2.B ■ B' csc(f2o^^)] j • 

(2.38) 

Once one takes care of a small technical issue, using this propagator in (2.37) gives the 
long-known result (1.5a). 

The technical issue has to do with appropriately handling a At —)■ 0 divergence in the 
integration in (2.37), which we will defer to section VIIA where a thorough discussion will 
teach us how to handle more substantial issues in the calculation of double splitting. 


III. REDUCING N PARTICLE TO N-2 PARTICLE QUANTUM MECHANICS 
A. Reduction from 2-dimensional point of view 

We now give a general argument about reducing N particle quantum mechanics to ef¬ 
fectively N—2 particles in the context of splitting rate calculations. First, we will identify 
a special sub-space of the Hilbert space that the A-particle Hamiltonian "H of (1.4) leaves 
invariant, in the special case that 

N 

^mi = 0. (3.1) 

i=l 

Then we will show that this is the relevant sub-space for our problem. 

As already noted, transverse translation invariance means that total transverse momen¬ 
tum is conserved, 

[W.E.PV). (3-2) 

SO that we may focus in particular on the subspace 

E,P. = 0 (3.3) 

relevant to our problem. We will show that one may simultaneously impose the constraint 
that 

= 0, (3.4) 

and that this also characterizes the relevant sector for our problem. From the form (1.4) of 

n, 

-i [n , Yi^ibi] = YiPi: ( 3 - 5 ) 

which vanishes in the center-of-mass frame (3.3). In ordinary quantum mechanics, the 
relation (3.5) just says that total mass times the time derivative of center-of-mass position 
equals total momentum. If the total momentum vanishes, as in (3.3), then [H, Yi "^ibi] = 0. 
However, this does not normally imply that the operator Yi^ibi is constant when in the 
center-of-mass frame because 


-i \Yi^ibi , EjPi] = 


(3.6) 
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bi2,Xi2 - b2,X2 

bii,Xii - bi,xi 

FIG. 13: The notation used in (3.10) for the state of the system just before and just after the first 
splitting in fig. 12. 


does not normally vanish, and so we may not simultaneously choose rriibi constant and 
J2iPi — 0- (3-1), however, (3.6) does vanish, and so both (3.3) and (3.4) may be 

imposed simultaneously. Time evolution with 'H will keep us in this sector if we start there. 

Do we start there? For the sake of concreteness, begin by considering the case of single 
splitting, as in £g. 12 and eq. (2.34). First consider the state of the system before any 
splittings. We start with |pi)(pi|, which in our language we would call a “2-particle” initial 
state 

|Pi,P2) = |-Pi,Pi) (3.7) 

in ]HI(i) ® IHI(i), with ( 0 : 11 , 0 : 12 ) = (—1,1). Because of the symmetry of the problem, we 
may formally average over the symmetry transformations (2.24) without changing the hnal 
answer. This replaces |pi)(pi| by something proportional to 



(3.8) 


which is 

[ \-PhPi) = f \buh) (3.9) 

Jpi Jbi 

in the language of ]H[(i) (g) EI(i). From the left-hand side of (3.9), note that we start in the 
sector J2iPi = Pi +P 2 = 0; from the right-hand side, we simultaneously start in the sector 
Y,i Xibi = Xiibi -F Xi2b2 = -bi + b2 = 0. 

Now we need to ask whether having o:iihii -|- Xi 2 bi 2 = 0 before the splitting at f = G in 
£g. 12 means that we will have Xibi +X 2 b 2 Toiaha = 0 just after the splitting. At t = G, the 
conjugate particle (particle 1) is just a spectator, and so (hi,a:i) = (hii,xii). (See £g. 13 for 
reference.) As to the remaining particles: the splitting vertex is local, and so b 2 = b^ = b^. 
Conservation of longitudinal momentum gives 2:2 + 0:3 = Xi 2 - Putting this together gives 


Xiihii -h a:i 2 hi 2 = Xibi -h X2b2 + Xsb^, (3.10) 

and so the condition of Yli^ibi = 0 is preserved by splitting. This argument works for any 
splitting with any number of particles, and so the condition is preserved by all the splittings 
in the double-splitting diagrams of figs. 5 and 6 as well. 

For the 3-particle case, the constraint Xihi 4 - 0 : 2^2 + 2 : 3 hs = 0 (along with xi -I-X 2 + X 3 = 0) 
explains the string of equalities asserted previously in (2.29), 

^ b2-b3 ^ ^ 

(xi + X2) (X2 + X3) (X3 -b Xi) 

Applying these constraints to the 3-particle Hamiltonian 'H of (2.11a) allows us to reduce 
l-i to the effective 1-particle Hamiltonian (2.31), as described previously. 


24 





B. Reduction from 3-dimensional point of view 


The constraint on Kiay also be understood in 3-dimensional language as a con¬ 

sequence of invariance under rotations that change the direction of the axis by a tiny 
amount, consistent with our large-p^ approximations. The conserved angular momentum 
is^o 

= ( 3 . 12 ) 

i 

where the superscript “(3)” indicates 3-vectors and S is the spin contribution to the angular 
momentum from the helicities of the particles. (In terms of the angular momentum J in 
the amplitude and J in the conjugate amplitude, S = J — J.) In the high-energy, nearly 
collinear limit, (i) the positions of the particles are all ~ t, and (ii) S will be in the ±z 
direction. Then the transverse piece JTI of (3.12) is given by 

- Y^ibiXiE - t YiPi^ (3-13) 

where J7’j_ = is the 2-dimensional dual of JT±- Since YiPi — 0) conservation of JT± 

then gives conservation of Yi^ibi (in the large-p^ limit). Since rotational invariance is a 
property of the full Hamiltonian (including the 5H for splitting), this gives conservation of 
YiXibi over the entire process. 


IV. DOUBLE SPLITTING USING JUST N-2 “PARTICLES” 

We now turn to writing expressions for contributions to double splitting, focusing hrst 
on the xyyx interference of £g. 11. Before we start, we should clarify something about the 
number of particles involved at various times in this diagram. 


A. N=5 particle evolution unnecessary 


The right side of fig. 11 is drawn in a way that indicates only three particles (two in 
the amplitude and one in the conjugate amplitude) are involved in the final time interval, 
ty < t < tx- However, based on the left-hand side of the figure, one might think there 
are five particles in that time region: all three final daughters {x, y, and 1—x—y) in the 
amplitude, and two particles (a final daughter y and an intermediate daughter/parent 1—y) 
in the conjugate amplitude. That is, the terms corresponding to the evolution during the 
time period ty < t < would be of the form 



'ffcfKf, txlpkK, ty) (PxKf, tx\pK, ty)' 


(4.1) 


For our purposes, the medium does not need to be invariant under large rotations. It is enough for it 
to be invariant over transverse distances probed during the formation time. For example, A6 ~ 1/Q± ~ 
foi- thick media, where w is the smallest energy of a daughter. 
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FIG. 14: Labeling in (4.1) for momenta at beginning and end of time interval tj<t< 


with momenta labeled as in fig. 14 and where we have also included the integration over the 
hnal-state momentum Kf of the y daughter. We can write this in 5-particle notation as 

/ (4.2) 

J K.f 

The 5-particle problem could be reduced to an effective 3-particle problem using the con¬ 
straints of the last section. 

Now recall that in the interval between splittings, the high-energy particles do not di¬ 
rectly interact with each other but instead simply propagate in the background fields of the 
medium, which are later averaged with ((•••)). So, when computing (pffcfKf, fy) in 

(4.1) for a particular background field configuration (before averaging), the three particles 
in the amplitude evolve independently with the 1-particle Hamiltonian if(i) appropriate to 
that particle. That is, 

= (pf,fx|P,%)(fcf,4|fc, %)(/«£, (4.3) 

and similarly for the two particles in the conjugate amplitude. We may therefore rewrite 
(4.1) as 

(Px,4|P,%)* y (Kf,fx|K,%) (Kf,tx|K,fy)*y (4.4) 

Inside the double angle brackets (before averaging), the 1-particle evolution of the y daughter 
is unitary, and so we may use the sum over the final states of the y daughter to simplify the 
above expression to 

<|{pffcf,fx|pfc,fy) (Px,fx|p,fy)*^^ {2Try6^‘^\K- k). (4.5) 

What is left is what we would call (medium averaged) 3-particle evolution (the x and 1—x—y 
daughters in the amplitude, and the intermediate 1—y daughter/parent in the conjugated 
amplitude), which can be reduced to effective 1-particle evolution using the method of the 
previous section. The lesson is that the sum over final states of a final-state daughter may be 
performed as soon as it has been emitted from both the amplitude and conjugate amplitude, 
and one need not worry about its interactions after that.^^ 
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There is a hidden assumption in this argument, which is that when we take the statistical average ((• • •)) 
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B. First expression for xyyx interference 


As mentioned before, the diagrams in fig. 6 may therefore be treated as 3-particle evolu¬ 
tion followed by 4-particle evolution followed by 3-particle evolution. After reducing each N 
particle problem to an (iV—2) particle problem, we may reinterpret each diagram as effec¬ 
tively 1-particle evolution followed by 2-particle evolution followed by 1-particle evolution. 
We will work in h-space with the variables 


B 




bj - bj 

Xi + Xj ’ 


(4.6) 


of which N—2 are independent for N particle evolution, once we impose the constraint 

Xibi = 0 . 

For now, we focus attention on the xyyx interference depicted by fig. 11 and the first 
diagrams in figs. 5 and 6. We will see later that the other interference diagrams can be 
related to the xyyx calculation, and so most of our task will be the evaluation of xyyx. For 
the iV=4 particle evolution, we label the Xi as shown in £g. 15, so that 


(xi, X 2 , £ 3 , £ 4 ) = (-1, y, ^-x-y, x). (4.7) 

(Our numbering convention is not quite arbitrary: avoiding having x and y cyclically se¬ 
quential in this list will turn out to be convenient for discussing the large-W limit.) For the 
sake of avoiding confusion, we find it convenient to put hats atop the Xi when specifically 
referring to the Xi that characterize the iV=4 evolution of xyyx, and we will often use these 
variables in other time intervals as well. For example, the three longitudinal momentum 
fractions associated with the initial 3-particle evolution for < t < ty are 

{xi,X2, X3) = (xi, -X1-X4, £4) = (xi,X2+X3, £4) = (-1, 1 -X, x), (4.8) 

and those associated with the final 3-particle evolution for ty < t < are 

(xi,X2,X3) = (-£3-£4,£3,£4) = (£ 1 +^ 2 , ^3, £ 4 ) = (-( 1 -y), 1 -x-y, x) . (4.9) 

For the sake of later relating results for xyyx to other interference contributions, it will be 
useful to keep many things in terms of the £i instead of writing formulas directly in terms 
of X and y. 

In order to help emphasize which variables apply to 4-particle evolution and which to 
3-particle evolution, we will refer to the Bij of (4.6) by the symbol in the 4-particle 
case. For the splitting matrix elements at the intermediate time ty, it will turn out to be 
convenient to choose our two variables describing the 4-particle (effectively 2-particle) system 


of the background fields we may ignore correlations of interactions just before the splitting time ty with 
interactions just after ty, within one correlation length of the plasma. This sort of assumption is implicit 
in many aspects of this work and standard treatments of the LPM effect, and its justification is that the 
effects of the background over a single correlation length about t = ty should be parametrically small 
compared to the effects over the entire formation time fform oc This is similar to the justification of 

footnote 5 for using an instantaneous potential F in "H, and the same sort of argument is implicit for our 
treatment of all the splittings in this paper. 
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—Xi—X4:,h x^,hz 



FIG. 15: Our labeling conventions for longitudinal momenta Xi and helicities hi for the xyyx 
interference diagram. 


to be (i) the Cij involving the two spectators to the splitting, and (ii) the Cij involving the 
two particles involved in the splitting. The same will be true at the other intermediate time 
ty. So, referring to £g. 15, the convenient 4-particle variables will turn out to be (C 41 , C 23 ) 
at ty but (C 34 , C 12 ) at ty. With these conventions, the xyyx interference contribution to 
double splitting has the form 


dl 

dx dy 

X {Cl,Cl,,ty\Cl,,Cl,,ty) 

X {B\ty\B^,U) {B^\-t 6 H\), 



(4.10) 


where all intermediate transverse momenta are implicitly integrated over—i.e. there is inte¬ 
gration over B^, B^, CI 4 , CI 2 , Ch, C 23 , B^, and B^ above. The superscripts are just labels 
for the different intermediate states (and are written as superscripts rather than subscripts 
for the sake of compactness of names like C^^). The initial N=2 (effectively Al—2=0) state 
is represented by |) and will be discussed later, along with the hnal state (|. The specihc 
particles involved in each 6 H splitting above should be inferred from the xyyx diagram in 
£g. 15; we have written dH above to indicate splittings that are in the conjugate amplitude. 
For comparison, the single splitting result, corresponding to £g. 10 and its conjugate, has 
the form 


dl 

— = 2 Re 
dx 


dx 


with 


dl 

dx 


E 


f J2^\t6H\B^) {B^U\B^,U) {B^\-zSH\) 

pol 


(4.11) 


(4.12) 


in this notation. 

We shall discuss the SH matrix elements in sections IV D and IV E. They have the form 


{B\SH\) = -^gr^^,k■V6^^\B) 
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(4.13) 













for parent i to split into daughters j and k. Here the two-dimensional vector 'Ti^jk 
form 


r 


i^jk 


n~'co\or'^ 


^jk 


2^3/2 


has the 

( 4 . 14 ) 


where is the color factor associated with the vertex [e.g. {T^)bc = —for gluon 

splitting g ^ gg with color indices a —)■ 6 c]. P is related to square roots of spin-dependent 
DGLAP splitting functions and depends on the Xi and the helicities of the parent and 
daughter particles in i —)■ jk (which also determine the helicity of the vector 'P). More 
details will be given later. 

The corresponding matrix elements for 5H linking the 3-particle to 4-particle sector will 
be reducible to the previous case (4.13). For example, if the daughters are particles 2 and 
3, we will hnd that 


= -igr^^23 ■ V(5(2)(C23) |x 4 + Xil-' - B). (4.15) 


Diagrammatic rules which summarize all the splitting matrix elements needed are given in 

£g. 16. 

Let’s now apply these rules to the xyyx interference depicted in fig. 15 for the case where 
X and y are bremsstrahlung gluons. In the high-energy limit, the interactions of high-energy 
particles with soft background fields preserve the helicities of the high-energy particles, 
and so helicities of individual high-energy particles are preserved by the propagators in 
fig. 15. The color charges of individual particles are not preserved, and so propagators like 
ty\C\^, CI 3 , ty) in (4.10) generically depend on color. As we’ll discuss, the large-A^, 
limit that we will take in this paper simplifies the color dynamics so that we may ignore this 
complication. Applying the rules of fig. 16 to fig. 15, naively contracting color indices, and 
using the 6 -functions to perform as many integrals as possible, then yields (see Appendix A 
for more detail) 


dl 

dx dy 


xyyx 




-1 


E 




B^ = 0 


X P-ft„ft_ft,(-£ 3 ,£ 3 +£ 4 ,-£ 4 ) ■ V,ty) 

X {V_-h^hu-hy{Xi+X2,-Xi,-X 2 ) ■ '^cl^){'P-h,h^,hy{-X2-Xz,X:i,X2) ' 

[f^y r^y +_\c'y c'y + \ 

0J2=U=0'23; V34=if^; 

(4.16) 


lc[2=o=c|3; cl^=By-, Cl,=By 
X 'P-hi,h,hA^l,-Xl-X4,Xi) ■VB-{B^,ty\B^,t^) 


B^=0 


where dn is the dimension of color representation R. For the case considered in this paper 
of an initial gluon, the color factor is 


d-ZiriUTiUTi) 


1(^2 _ 1 Ar2 
2^A — 2-^^c ■ 


( 4 . 17 ) 


C. Large-Afc limit 

When discussing the single splitting rate in section II, we were able to determine the dot 
products Ti ■ Tj of color charges in the effective potential (2.18) in terms of fixed Casimirs 
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) ®A: 




FIG. 16: The diagrammatic rules for splittings linking (via either —i6H or +i5H) the state |) to 
\B) (top rule) or \B) to 1 ( 734 , Cu) or permutation thereof (bottom rule). Buv = {bu — bv)/ixu+Xy) 
and may refer, in different contexts, to ± the 3-particle B, or one of the 4-particle Cuv, or to some 
mixture. However, note that Bji = Bkj = Bik in the top rule, which can be used to always 
write expressions in terms of 3-particle B and/or 4-particle Cij’s. The blue arrows on the particle 
line indicate color flow of color representation R. (In the case of R=A, appropriate to g ^ gg 
splitting, the direction of the color flow does not matter.) bi, ai, xi, and hi indicate the transverse 
position, color index, longitudinal momentum, and helicity of each particle. The black arrows give 
the convention for the flow of xi and hi in the statement of the rule, and these values should be 
negated if they are instead defined by flow in the opposite direction. In the bottom rule, color and 
helicity indices and their contractions are not explicitly shown for the spectators because they are 
trivially contracted. Conservation of longitudinal momentum means Xi + Xj + Xk = 0 (top) and 
additionally Xm = Xr and Xn = Xg (bottom). 


Ci, and so we avoided having to deal with color dynamics. This simplihcation is no longer 
possible, in general, when analyzing the N=4 particle evolntion that appears in donble- 
scattering: Ti -|- T 2 -|- T 3 -|- T 4 = 0 is not snfiicient to determine all the Ti ■ Tj for 4 particles. 
In this work, we sidestep this complication by focnsing on the large-iV/ limit. 

Becanse of the kinematics of onr problem (high energy particles interacting with the 
plasma via momentnm transfers that are soft compared to E), it has been convenient in onr 
analysis to time-order onr interference diagrams, as in £g. 6 . To discnss the large-limit, 
however, we wonld also like to draw onr diagrams as planar diagrams. In the case that all 
particles are glnons, the diagrams of fig. 6 can be drawn as time-ordered planar diagrams 
by drawing them on a cylinder (since a cylinder can be mapped into a plane), as shown in 
fig. 17a for the xyyx interference. 

Fig. 17b shows the color flow of £g. 17a decorated by an explicit example of backgronnd- 
field interactions and correlations in the large-iV/ limit. The reqnirement of planarity means 
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time 


time 




A B 



xyyx 

(a) 



A __ B 

xyyx 

(b) 


FIG. 17: (a) The xyyx interference diagram drawn on a cylinder. The top edge AB of the shaded 
region is to be identified with the bottom edge AB. (b) The same showing color flow for an explicit 
example of medium background field correlations (black) that gives a planar diagram (and so 
leading-order in 1/Nc). Bonus craft project: Cut out these figures with children and glue into 
cylinders. Collect and trade with friends! 


that each fundamental-charge high-energy color line shown in £g. 17b only interacts (via 
medium correlations) with one other, and those interactions are independent of the in¬ 
teractions of other color lines. That is, the general A^-particle potential in 'H among our 
adjoint-charge particles reduces to 


V{bi, - ■ - bNit) 


Vp{b2—bi; f) -f-Vp(h3—^2; t) + - ■ ■ + Vp{bN—bN_i;t) + Vp{bi—bN; t), (4.18) 


where Vp is the large-A), potential (in the language of H) between a fundamental charge and 
anti-fundamental charge particle, and hi, 62 , b^ label the high-energy particles cyclically 
as one goes around the cylinder (e.g. from bottom to top in fig. 17a). In the harmonic 
oscillator approximation, this becomes 


—^ (^2 ~ biY + (^3 — ^ 2 )^ + • ■ ■ + (hjv — b]\f-iY -|- (hi — h 
which (in the large-A), limit) is the same as 

hb2 - biY + (hg - hs)^ -b ■ ■ 


'N 


{bjsf — hjv-i) + (hi — h 


'N 


(4.19) 

(4.20) 


Note that for N=3 this agrees with the all-gluon case of the more general formula (2.21), 
which did not depend on large Nc- 


D. Reduction of states from N to N—2 particles 

We now turn to an explicit construction of the N—2 particle states such as |C 34 ,Ci 2 ) 
appearing earlier. The utility of an explicit construction is that it will allow us to hnd the 
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correct normalization factors associated with 6H matrix elements, like the \xm + in the 
bottom-right of £g. 16. 

For = 0, we may project any b-space state onto the subspace with = 0 

and YhiPi = 0 by (i) imposing Yhi^ibi = 0 in b-space, and (ii) averaging over all b-space 
translations in order to project onto YliPi ~ b- Treating all N particles symmetrically, this 
procedure defines 


\{B,^})=Af[ |bi + Ab,---b^ + Ab) 

JAb 


Y,Xibi=0 


(4.21) 


where represents the (not independent) set of all of the values of Bij = {bi — bj )/(xj -|- 

Xj). M is an Xj-independent normalization factor described in appendix B. In the appendix, 
we show that this particle-symmetric dehnition of the projected states normalizes those 
states as 


({CylKC',}) = (43 + (N=i), (4,22a) 

({B«}|{B',}>=iPI(Bi2-B;2) (A'=3), (4.22b) 

(|)=a; (1V=2), (4.22c) 


where we have used our convention of referring to the 4-particle variables as Cij and Xi. One 
may use any permutation of indices on the right-hand sides above. Note that (£3 + £ 4 )^ = 
{xi -f- £ 2 )^ in the A=4 case and x\ = x\ in the N=2 case. 

The notation IC 34 , C 12 ) that we used earlier is related to the above symmetrically-dehned 
states by 

10*34, C 12 ) = |a;3-l-a;4| 1(0^}), (4.23) 

so that the A=4 states are more conveniently normalized as 

{C 31 , C„\C' 3 „ C[,) = (4.24) 

The N=3 state |.B) used earlier is simply 

|B> = |{B«}>. (4.25) 

The normalization (4.22c) of |) may seem a bit strange, but |) is what we want for the 
hnal state. At the end, after the emission of both the x and y bremsstrahlung gluons, we 
want to sum over all possible final states of the emitter, 

f \Pi){Pf \ = f l^f)(bf| (4.26) 

^ Pf J b{ 

which in our A-particle notation is 

lend) = f |bf, bf). (4.27) 

Jbi 

By a shift of integration variable, this is with |) defined by (4.21) as 


\)^M 



bi-^Ab, b2-FAb) 


x\hi+X2b2=0 


(4.28) 
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noting that Xihi + X 2 h 2 = 0 enforces bi = 62 , given that xi + X 2 = ^- A similar argument 
can be made about the initial state if we hrst use the symmetry of the problem to average 
over small changes in the directions of the initial momentum p\ as in (3.7-3.9). Details of 
how the normalization factor M cancels in the hnal result are given in appendix B. 

In the next sub-section, we will discuss in detail the form of the 5H matrix elements for 
individual particles. For splitting i ^ jk in the amplitude, this has the form 


{Pj,Pk\5H\pi) = gTi^jk ■ Pjk (4.29) 

with Pjk = XkPj — XjPk-^"^ We then need to project this result from N particle language 
to N—2 particle language. In appendix B, we show that (4.29) implies, in our constrained 
subspace, the results (4.13) and (4.15) quoted previously for the matrix elements {B\SH\) 
and {Cu,C 23 \SH\B). 


E. Splitting matrix elements 

The matrix elements for nearly-collinear bremsstrahlung (or pair creation) is text¬ 
book material, used, for instance, in discussions of DGLAP evolution. In our case, we 
will need individual helicity-dependent matrix elements (e.g. from ref. [24]) and not just 
helicity-averaged matrix elements. We also need a few conversions compared to the stan¬ 
dard literature. Our states \p) have been normalized with non-relativistic normalization 
{Pi\p'i) = (271)2^(2) because we have cast our problem into the language of non- 

relativistic quantum mechanics. Results in the literature are instead for states with rela¬ 
tivistic normalization {pi\p'i) = 2Ei{2TT)‘^6^'^\pi — p') [times 27r 6{piz — so the 

matrix elements differ by factors of \/ 2 A) = y/2\xi\E. Secondly, results in the literature are 
typically specialized to a choice of 2 : axis that is exactly aligned with the direction of the ini¬ 
tial particle: p* = 0 for p* —)■ PjPk- We will want to generalize to slightly different directions 
of the 2 : axis and so will express results in terms of the invariant Pjk = XkPj — XjPk instead 
of, say, pk- Finally, results in the literature are generally expressed with the convention that 
E is the energy of the parent particle. In our problem, the parent of the second splitting is 
a daughter of the hrst splitting, and so we will write parent energies more generally as \xi\E 
instead of E. Doing this will also make manifest in our formulas the symmetry between the 
three particles (parent and daughters) participating in a single splitting process. 

We are going to just present here the formulas for the matrix elements. Readers interested 
in details about the connection to standard formulas in the literature should see appendix 
C. The matrix elements are given by (4.29) with^^ 

^color-^ 

- 7 - _ i^jk t^jk , , 


22 Note that Pjk = Ph = Pij if we were to pick the convention that all xi flow away from the vertex, so 
that Xi is negative, and were to similarly negate Pi. 

2 ^ In this paper, we have generally swept the longitudinal momentum conservation factors 2t: 5(piz — Piz) 
under the rug. See appendix A on (4.29-4.31) for a more detailed discussion. 

2 ^ See appendix A for a note on dimensional analysis. 
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and 


i^jk 


'(di) •) ^k') 


\XiXjXk\ 


C 


R 


where the Pi^jk are helicity-dependent DGLAP splitting functions and 

G(^) = G3. ib 


(4.31) 


(4.32) 


are ± helicity basis vectors in the transverse plane with sign determined by the conservation 
of the (^-component Jz of angular momentum in the nearly-collinear limit. This is easiest to 
state if we adopt the convention that each helicity flows away from the vertex. So, instead 

of thinking of the helicities in a splitting as, for example, H-)■ ++, we will think of them 

as (hj, hj, hk) = (—, +, +). With this convention, (4.31) is, more precisely. 


i^jk 


^{hi-\-hj-\-hk) / Phi,hj,hk i^i') 5 ^k) 


\XiXjXk\ 


c 


R 


(4.33) 


For g ^ gg with Xi —)■ XjXk, the DGLAP splitting functions (for \xj\, \xk\ > 0) are 


P-i 


= Gai 


Xi 


\XiXjXk\ 


= Gai 

p+^- = 0 , 


Xi 


\XiXjXk\ 


P+^_+ = Gai 


Xl 


\XiXjXk\ 


(4.34) 


corresponding to 



Xi 


XiXjXk\^l‘^ 

e(-)l 


Xi 


\XiXj 




XI 


\XiXjXk\^/‘^ 


e(-)l 


V+^ = V = 0 . 


(4.35a) 

(4.35b) 

(4.35c) 

(4.35d) 


The remaining cases are obtained by negating all helicities and the subscript of e(±). Note 
that the 7^’s do not depend on the signs of the xi. Note also that in the soft bremsstrahlung 
limit \xk\ \xi\, helicity is conserved among the harder particles (i.e. H—)■ ++ and H—)■ H— 
are unsuppressed but H—)■ —h is suppressed) and |P+^.++| ~ |P+^.+_|. These simplihca- 
tions account for the fact that a careful treatment of intermediate helicities was not generally 
necessary in earlier work [6-8] focused on the soft bremsstrahlung limit. (See Appendix F2 
for a caveat.) 


F. Helicity sums 

It is now time to do the helicity sum in our expression (4.16) for the xyyx interference. 
The only things that depend on helicities are the P factors. These factors will be the same 
for all of the diagrams in the first row of hgs. 5 and 6, and in this case it will be more 
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convenient to write in terms of x and y explicitly rather than in terms of the 4-particle XiS. 
Specifically, we need 

/lx 5 /^y 1^'z 

X (4-36) 

where P” are the Cartesian components of 7^. If one prefers, this can be written equivalently 
as 

^K^h.,hS^-y ^ ^ i-y,!/)] 

hx:hy,hz h 

X [J2'^h^h,,hy{^-^ ^ ^-x-y,y) l-x,x) , (4.37) 

h 

which (reading right to left) shows more clearly how the successive splittings tie together in 
both the amplitude and conjugate amplitude. 

Because of parity/reflection symmetry in the transverse plane, the rate (4.16) is inde¬ 
pendent of the initial helicity hi, and so we may replace (4.37) by its helicity average over 
hi. Then, because of rotation invariance of the transverse plane, the above helicity sum 
(averaged over hi) must have the form 

a{x, y) + f3{x, y) 5^^5^^ + 7 ( 2 ;, y) ^- 5 -- (4.38) 

for some functions a, (3, and 7 . Using (4.35) in (4.37), we hnd 



y^{l—x)^{l—y)^{l—x—y) x^{l—x)^{l—y)^{l—x—y) 


1—x 1—y 

x^y^{l—y)^{l—x—y) x^y^{l—x)^{l—x—y) 


_I_I_f_ 

y^{l—x){l—y){l—x—y)^ x^{l—x){l—y){l—x—y)^ 

l—x—y 1 

^ x^y^{l—x){l—y) x^y^{l—x){l—y){l—x—y)^ 

1 —a; 1—y 

xy{l—y)^{l—x—y)^ ^ xy{l—x)^{l—x—y)^ 

l—x—y 1 

xy{l—x)^{l—y)^ xy{l—x)^{l—y)^{l—x—y)^ 

for the case studied in this paper, where all high-energy particles are gluons. Substituting 


(4.39) 
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(4.38) for (4.36) in the xyyx expression (4.16) gives 


dl 


( 

dx dy 

- SE^ 

xyyx 1 

ail -F £4 £3 + £4 

J 


' By,By 




B^=0 


X ^cfycyJCL CL ty\Cl, Cl,, t,) 


Cf,=0=C,^,; Cl.=By-, Cl,=By 




B^=0 


(4.40) 


Note that a, (3, and 7 are all symmetric under x -H- y. 

This is as far as we can go without either doing numerics or making additional approxi¬ 
mations, to which we now turn. 


V. xyyx IN HARMONIC AND THICK-MEDIA APPROXIMATIONS 

We can make further analytic progress by (i) restricting attention to the harmonic ap¬ 
proximation, and (ii) assuming that the medium is sufficiently homogeneous that q does 
not vary signihcantly over distances of order the formation length. In this section, we will 
continue to focus for now on the xyyx interference. We will tackle the iV=3 particle ends of 
the evolution in fig. 15 and integrate over the hrst and last splitting times, C and Then 
we tackle the iV=4 particle piece of the evolution by solving the harmonic oscillator problem 
for the propagator in the large-W limit, and then we do the B integrals remaining in (4.40). 
This will leave the result in the form of a single 1-dimensional integral over At = ty — ty. 


A. Integrating over first and last times C and 

During the 3-particle phases of the evolution of the xyyx interference of fig. 15, our "H is 
given by (2.33) and (2.27). For the initial 3-particle evolution corresponding to < t < ty, 
we’ll write this as 

TD'2. 

(5.1) 
(5.2a) 

(5.2b) 




with 


Mi = XiX 4 ,{xi+X 4 )E = x{l — x)E 
and (when all particles are gluons) 


—— 1 

2E Vail £4 a;i-|-a; 4 / 




2E 




X 1—x 


For the final 3-particle evolution corresponding to E < t < t^. 


n = 


2 Mf 




with 


= (ai 3 +a; 4 )a; 4 a; 3 i 7 = (1 — y)x{l — x — y)E 


(5.3) 

(5.4a) 
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and 


a = w-^ 

2E 


4- + I - / 

+ - + -I -W-^ 


11 1 

- 1 -^ ^1 - 

1—y X 1—x—y 


X3 + X4 X4 X3 

(Throughout this paper, one should understand \/S to refer to the root 
The propagator (2.38) for a harmonic oscillator is 

Mhl csc(r2t) 


(5.4b) 


{B,t\B',0) = 


2Tri 


exp (I Mhl {B"^ + B'^) cot (hit) — 2B ■ B' csc(r2t)] 


The factor corresponding to the initial 3-particle evolution in (4.40) is of the form 

csc^(r 2 (t — t')) 

VB'{B,tB',t') = -^^ 

B'=o 2n 

Integrating over the initial time t' < t above gives 


B exp UMVL cot {VL{t - t'))B‘^‘ 


dt' VB'{B,t\B',t') 


B'=0 


iMB 

ttB^ 


exp 


cot {n{t - t')) 


t'=t 


t'=—oo 


(5.5) 


(5.6) 


(5.7) 


The term with t' = t has a divergent exponent, corresponding to inhnitely oscillatory be¬ 
havior in B. Naively, one might expect to be able to drop it, since an inhnitely oscillatory 
function vanishes if later integrated against a smooth function. We will drop it here but will 
later have to return and be more careful about what happens when times become coincident. 
[In particular, if a third time in the double splitting problem becomes coincident with t and 
t' above, then (5.7) is not integrated against a smooth function of Bt\ 

In our problem, hi has a negative imaginary part, so that cot(r2cx)) = i, and so the 
t' = —oo term in (5.7) leaves us with 

r*t 


f dt’VB'{B,t\B',t’) 

J —OO 


B'=0 


iMB 

7iB‘^ 


exp{-lMQB^). 


(5.8) 


For the sake of evaluating other interference diagrams in the future, it will be useful to 
generalize to include cases where hi may have a positive imaginary part. See appendix A for 
an argument that the more general result may be cast in the convenient form 

ct 


[ dt'VB'{B,t\B',t') 

J —OO 


B'=0 


iMB 

nB"^ 


exp 


■l\M\nB^). 


(5.9a) 


A similar analysis, relevant to the hnal 3-particle evolution, gives (with the same caveats) 


r+oo 

/ dtVB{B,t\B',t') 

Jt' 


B=0 


iMB' 

' 7r5'2 


exp(-i|M|f25'2). 


(5.9b) 


Using (5.9) on the expression (4.40) for the rate gives 
dx dy 


xyyx 


X 


•OD 

d{At) 


\xi + a;4||a;3 -7 £4 

f Bt Bl 

'By,By {B^Y {B^Y 


exp(-i|Mf|a(5")2 - \\M,\n,{B^f) 


^ <^" 2 , AtICli, Cf 3 , 0 ) 


C'b=o=a3; C|4=By Cl.=By 


, (5.10) 
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where At = ty — ty, and we have used time translation invariance (in our thick-media 
approximation) to convert the expression for the probability I into an expression for the 
rate F. 


B. 4-particle normal modes and frequencies in large Nc 

We now turn to the iV=4 particle propagator (CI 4 , CI 2 , At\Cl„C^23,0). 


1. Reduced Lagrangian 


Consider the 4-particle Lagrangian 

L = \mih\ + \ni2h\ + — V (5-11) 

in the translation invariant case that V depends only on differences = hj — bj of the 
coordinates. In the particular case where ~ 0; want to impose the constraint 

rriibi = 0 on L above. The result can be written in the form 


L = 


-h lm2bM + -V{bi^, ^24, ^ 34 ) 


mibi=0 


(5.12) 


as can be seen by expanding the factors of (hj — 64 )^ and using the constraint and Yhi = 0 
to recover (5.11). The constraint may be rewritten in the form mihi 4 -|-m 2 h 24 +^^ 3 h 34 = 0, 
and so we may rewrite L in terms of just two variables as, for example, 


L = 


i mibl^ + i m2bl^ + | msbl^ - V(bu, ^24, ^ 34 ) 

T 


f> 24 =—(mi 614 +7713 634)/m2 


^14 


^34 


M, 


( 2 ) 


hi4 


'34 


— V(bi4, b24, b 


'34 j 


f> 24 =—(mi 614+m3 634)/m2 


(5.13a) 


with 


M(2) = 


_ mi 


m 3 


mim3\ ^ f ^ ^ 


m2 \m3m1 m3 


1 

7714 


— + — 
7714 


-1 


(5.13b) 


We could now solve this 2-particle problem. However, it will be convenient to hrst change 
to one of the sets of variables that we actually use in our rate expression (5.10). We will 
trade the hu and h34 used above for C12 = bi2/{xi + X2) and C34 = b34/{x3 + X4). Using 


25 Pqj. readers wondering how we could get from the 4-particle description to the 2-particle description 
without ever explicitly imposing the condition P 1 -I-P 2 +P 3 +P 4 = 0 for our subspace: Whenever we can 
consistently project to ^^rriibi = 0 for all times (as demonstrated by the earlier arguments of section 
III A), then = 0 follows automatically from (3.5). 
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= 0 and Xibi = 0 , one can show that 


hi 2 — ( 2^1 + X 2 )Ci 2 ) (5.14a) 

hi3 = X 2 C 12 — (5.14b) 

hi4 = X 2 C 12 + x^Cm, (5.14c) 

^23 = —a^iC '12 — a; 4 G 34 , (5.14d) 

^24 = —X 1 C 12 + 0 : 3034 , (5.14e) 

^34 = ( 3:3 + 3 : 4 ) 034 . (5.14f) 

Then (5.13) becomes 

^ ^ (tu) “ ^ 


with 

A nice feature of these variables is that the kinetic term is diagonal, unlike in (5.13). Note 
that, for the Xi of (4.7) relevant to the xyyx interference, 9H above is positive dehnite. 

In harmonic approximation, (5.15) will reduce to the form 

‘=Ka)’-(a)-Ka)'*(s;)- 


with ^ depending on specihc details that will be given in a moment. The squares r2| 
of the normal mode frequencies will be given by the eigenvalues of 011“^/^. The 

corresponding normal modes (O34, O 42 ) will be orthogonal with respect to Olt, and we will 
normalize them as 

We will now turn to specihcs by deriving normal mode results in the large-W approxi¬ 
mation when all particles are gluons. But we will then package our results so that we can 
continue to simplify our rate expression (5.10) in a way that is independent of those details, 
so that our results can easily be adapted at a future date to other situations. 


2 . Specifics for large-Nc gluons 
Combining the large-W harmonic potential 

-^(642 + + ^ 41 ) (5.19) 

of (4.20) for gluons with (5.14) gives 

Q f 2 {xl + X3X4 +xl) X2X3 + X1X4, \ . . 

4 \ X 2 X 3 + X 1 X 4 2{xl + X 1 X 2 + xDj ' \ ■ ) 
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in (5.17). Note that, since this is a harmonic oscillator problem, the two transverse dimen¬ 
sions decouple, and so this may be treated as a problem in one transverse dimension rather 
than two. The result for the normal mode frequencies should have the same symmetry 
with respect to exchanging particle labels as (5.19) has, and (mostly for aesthetic reasons) 
we use = 0 to algebraically manipulate the result into a form where this is manifest. 

We hnd 


= 


' AE 


1 1 1 1 , 

-1-1-1-i vA ) 

Xi X2 X3 X 4 J 


1/2 


where 
A = 


XI 


Xn 


x\ xl X1X2 


1 


1 


> 0 


X2X3 X3X4 X4X1 X1X3 X2X4 

A more compact (but slightly less manifestly symmetric) expression is 

^ _ 1 1 1 1 (a;3-l-a;4)^ + {xiEx^f 

^ ~ „2 


X 


Xo 


Xn 


Xa 


X1X2X3X4 


With normalization (5.18), the normal modes (C^, C 12 ) and (C 34 ,Ci 2 ) are^® 

1 


C^34 = 


*-"12 — 


X2 


X3 + X4 
X4 


X1X3 

2N±E 


1 1 
- 1 -h 

X3 Xi X4 X3X2 


X1X3 


Xi+X2y 2N±E 


1 1 1 

Xi X3 X2 X1X4 




with 


N± = -a;ia;2a;3a;4(a;i - 1 - a; 3 )A ± (a;ia;4 -t- a;2a;3)(a;ia;2 + X3X4)\fK. 


(5.21) 


(5.22) 


(5.23) 

(5.24a) 
(5.24b) 

(5.25) 

For the Xi of (4.7) relevant to the xyyx interference, the xiX3/2N±E under the square root 
in (5.24) is positive. 


3. General form 


Whatever the specihc formulas for the normal modes are in a given problem of interest, 
we can write any (C 34 , C 12 ) as a superposition 



(5.26) 


of normal modes with superposition coefficients A^. and A_. We will write this in matrix 
form as 


C34\ 

C 12 ) 


— tty 



(5.27) 


The right-hand sides of (5.24) are in fact symmetric under 12,34 —>■ 21,43 (which takes N± —>■ —Nip), 
but we have not found an appealing way (i.e. other than brute force symmetrization) to write them that 
makes this symmetry manifest. [Also, one might think that 12,34 —>■ 21,43 should take (C'^,C'j^) —>■ 
iCt3,Ct) = (-C3±,-Cp±) instead of to (TC'^,+(7)^), but the overall sign of a normal mode is an 
arbitrary normalization choice, reflected by the ambiguity of the square root in (5.24).] 
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with 


Oy — 


= / ci ci 
C+ C.! 


(5.28) 


We use the subscript y because (C34, C12) are the variables we want to use at ty in (5.10). 

Because of the normalization (5.18) of our modes, the Lagrangian (5.17) in terms of A± 
is simply 

^ = E [5 'it - 5 . (5.29) 

± 

The corresponding propagator is 


(A+,A_,t|A'+,A'_,0) = 

csc(r2±t) 


n 


2m 


exp(i[i(A^ + A'|)n± cot{il±t) — A± ■ A'_j_r2± csc(r2±t)] 


. (5.30) 


At t = ty in (5.10), the variables we want to use are ( 041 , 023 ) instead of ( 034 , 012 ). 
However, it is easy to convert: (5.14) and Cij = bij/{xi+Xj) give 


O 41 

O 23 


_ -X 3 f O 34 

{Xi + Xi) V ^4 XiJ \^Oi2 


and so 


with 


<=^41 , _ 

C.’" ' y 


'23 


A 3 

A 


tty — 


-X 3 -X 2 


{X1+X4) V ^4 XiJ y 
The changes of variables (5.27) and (5.32) give^^ 

{Al,Al,At\Al,Al,0) 


(0f4,0[2,At|0Ii, 02^3,0) = 


det Ov 11 det 


From (5.28) and the normalization (5.18), ajTlay = 1, and so 

(which we will hnd useful later) and 

IdetOyl”^ = |detT)T|^A = \xiX2X3X4\^^‘^\x3+X4\E. 

By symmetry, or from (5.33), similarly 

IdetOyl”^ = 1x1X2X3X41 ^^"^\xi+X 4 \E, 


(5.31) 

(5.32) 

(5.33) 

(5.34) 

(5.35) 

(5.36) 

(5.37) 


so that 

(OI 4 , 0^2, At|OJi, 0^3,0) = \xiX2X3X4\\xi+X4\\x 3+X4\E‘^{A^_^_, AL, At\A\, AL,0). (5.38) 


See Appendix A for comments. 
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Combining with (5.27), (5.30), and (5.32), 


(CLCi2,At|CI„Cf3,0) = 


(27ri) {—XiX 2 X^X/^)\xi+Xi\\x^+Xi\E^VLj^VL_ csc(r2+At) csc(r2_At) 


X exp 


^yM Oy cot(fi At) a 


-1 (C^ 


y 


41 
/~^y 

*-"23 


I 

+ 2 
— i 


C'|4^ „-lT 

c'L 

r^y 
*-"41 
r^y 
*-"23 


Oy ficot(fiAt)a 


-1 


a ficsc(fiAt)ay 


-1 

y 


r^y 

*-"34 

r^y 

*-"12 

(~iy 

*-"34 

r^y 

*-"12 


(5.39) 


where 

f! ^ , (5.40) 

We have used (4.7) to rewrite \xiX 2 X 2 ,Xi\ above as —X 1 X 2 XSX 4 for reasons that are described 
in appendix E 1 b, having to do with relating the result for the xyyx interference to other 
interference contributions. 


C. Doing the B integrals 


We will now assemble our results and then carry out all remaining integrals except for 
the At integral. 

We want to use the propagator (5.39) in the expression (5.10) for the rate. It is convenient 
to hrst combine all of the exponential factors (before taking derivatives and setting various 
Cij to zero) and rewrite them as 


exp 




Ay Ay 
Ay Ay 



1 

2 





where^® 


I") 

V'-^23/ V^yy ^yyJ 


Ay Ay 
Ay Ay 

Ay Ay 
Ay Ay 


^yy ^yy 

Y - 7 - 
■' yy ^yy 


0^ -iay^’^ficot(fiAt)ay\ 
0^ _ ia-iT^cot(flAt) 
—ia~^^Qcsc{flAt) 



(5.41) 


(5.42a) 

(5.42b) 

(5.42c) 


Using (5.35), one may rewrite the factors of Oy ^ in (5.42) as Oy ^ = aJfSl, if desired. A similar relation 
between 0 “^ and Oy holds if one cyclically permutes (xi,X 2 ,X 3 ,a; 4 ) —>■ (x 2 ,X 3 ,X 4 ,xi) in the expression 
(5.16) for m. 
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With this notation, the rate (5.10) becomes 


dT 

\_dx dy 


xyyx 


ClajMM 

32n^E‘^ 




X / (i(At) r2+r2_ csc(r2+At) csc(r2_At) 

Jo 

r _ 

X [{YyBy - Y,yB%{YyB> - YyyB^U + ZyySnyy] 


X exp 


- \X,(BYf + X„B’ . B 


(5.43) 


The B integrals may be performed using 


Jo = / exp 

J By,By 


^X,{B’f - \X,(B^f + Xy,B^ . B 


47r^ 


h = 


(X,X, - X^,) ’ 

B^ ■ fiS 


exp 


'By,By {By)‘^{By)'^ 


jXyiByf - lx,(Byy + XyyBy ■ 


Xyy 




A = 


(By -By)^ 


'By,By 


(By)^(Byy 


■ exp 


-^Xy(Byy - \Xy{Byf + XyyBy ■ By 


27r^ 




In ( 1 - 


^y^y 


+ 


47r^ 


(XyX, - AT) ’ 


h = 


' By,By 


By By 


exp 


-^.XyiByf - \Xy(Byf + XyyBy ■ By 


AT^^Xyy 


Xy(XyXy - Xyy 


2 ^ ’ 


h = 


' By,By 


By By 

~Wf 


exp 


\Xy(Byf - \Xy(Byf + XyyBy ■ B 


AT^^Xyy 


Xy{XyXy - A ^ ) ’ 


(5.44a) 


(5.44b) 


(5.44c) 


(5.44d) 


(5.44e) 


giving 

dr 


dx dy 


xyyx 


CWM-M . . . . . 5 a 5 ^ ^ ... A 5 ... A 5 

:-— {~XiX 2 X 3 X 4 ) / d(At) r2+r2_ csc(r2+At) csc(r2_At) 

327r^ii/ Jq 


X 


yYy + aYyyYyy)lQ T (Cl T /3 T 2'^) Zyy 11 

+ [(o+7)yyyy + (/? + 7)yy5y,y]c - («+/3+ iW^y,!.^ + yy^yC)}. 

(5.45) 
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D. Small-time divergence 


The At integration in the xyyx result (5.45) has both a linear and a log UV divergence 
associated with At —)■ 0. Specihcally, if we expand the integrand in powers of At (see 
appendix D 1 for a little more detail), we hnd something of the form 


dV 

dx dy 


d{At) ■ 


stuff iC\al 


+ 


xyyx 

. 2 - -2 


(At)^ IhTT^At 


(rij sgn Mi + ftf sgn Mf) 


X xlx2xlx/^{xi -f X4y{x^ + x^y 


[oi -|- y') 


{a + '-))x 2 Xa 
{X i+X^yx^+XA) 


+ 0((At)°)l. (5.46) 


We have not bothered to explicitly show the coefficient “stuff” of l/(At)^ because that 
divergence is easy to dispense with. Here is one argument that is easy to make now. The 
coefficient “stuff” turns out not to depend on any of the frequencies H and so does not 
depend on q. It is therefore the contribution to (5.46) that comes from radiation in vacuum; 
it gives the result we would have gotten if we had set q to zero. However, the total rate for 
bremsstrahlung in vacuum must be zero by energy and momentum conservation, and so^® 


dr 


dr 


dr 

dx dy 

total 

dx dy 

total 

dx dy 


(5.47) 


total 


Thus no harm is done if we subtract the vacuum piece for each individual contribution that 
we compute and so compute 


dr 


dr 

dx dy 

xyyx 

dx dy 




xyyx 


(5.48) 


instead of [dV/dx dy]xyyx- This simply subtracts out the l/(At)^ piece of (5.46). (We will 
give an alternative argument later that the l/(At)^ divergences cancel without relying on a 
priori knowledge that the total vacuum contribution must vanish.) 

In contrast, the 1/At divergences in (5.46) do depend on the medium, but in a very specihc 
way: they are the sum of terms which depend on the medium (i.e. depend on q) either (i) 
only via Hi or (ii) only via Hf. In the first case, regarding the Q^/At terms, this means 
that only the contribution from the evolution before ty in £g. 15 depends on the medium, 
and so (regarding the same terms) the 3- and 4-particle evolution after ty is equivalent to 
vacuum evolution. One expects the medium to be irrelevant over very short times, and 
indeed we show explicitly in appendix D2 that the Q-JAt divergence arises specifically from 
the limit where ty, ty, and t^ approach each other simultaneously in fig. 15. Similarly, the 
Of/At divergence arises when tx, ty, and ty approach each other simultaneously. These two 
situations are depicted in fig. 18 for future reference. 


The vacuum rate is zero because we are analyzing the thick medium limit, which can be thought of as 
an infinite medium. If we had instead restricted splitting time integrals in (2.1) and (4.10) to start at 
some initial time 0, then there would be vacuum radiation associated with the initial appearance of the 
particle, related to the vacuum radiation associated with the hard process that created our initial particle 
(except lacking an ultraviolet cut-off). The usual analysis technique in this case would be to subtract off 
the vacuum result in order to isolate the part of dl/dx (or in our case dl/dx dy) due to medium effects. 
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M divergence 


-x-y-y -X 

divergence 


FIG. 18: A depiction of the (i) Hi/At and (ii) Of/At divergence of the xyyx interference as arising 
from the simultaneous approach of (i) ty, ty, tx versus (ii) tx, ty, %• These divergences correspond 
to the short-time approximation of vacuum evolution in the shaded regions. 


We will later see that the 1/At divergences cancel among different interference contribu¬ 
tions, and so we will be able to carry out the At integral numerically if we hrst combine 
results for different interference contributions before integrating. However, as will be dis¬ 
cussed later, we will nonetheless have to be careful to account for an additional contribution 
associated with the pole at At = 0. 

VI. OTHER CROSSED INTERFERENCE CONTRIBUTIONS 

We now turn to relating the other crossed interference diagrams in £g. 6 to the result just 
derived for the xyyx interference. Details are given in appendix E, and here we will instead 
loosely motivate and then describe the results. 


A. xyyx interference 

The second diagram of hgs. 5 and 6 describes what we call the xyyx interference con¬ 
tribution, which we depict with time labels in £g. 19. One difference with the xyyx case 
analyzed previously is that the Xi during the period of 4-particle evolution are different. In 
contrast to (4.7), we have 


(x[, X 2 , x'^, x'^) = {-(l-y),-y,l-x,x), (6.1) 

numbered (for the sake of large-W) in the cyclic order one would get if drawing the diagram 
on a cylinder, analogous to the xyyx case of £g. 17. Note that there are now two negative Xi 
because two of the high-energy particles are in the conjugate amplitude during the 4-particle 
evolution {ty <t < ty). We have chosen to indicate the 4-particle Xi for xyyx using primes, 
as in (6.1), and we will reserve the hats for the xyyx case, as in (4.7). The two are related 
by 


{x\,x'2,x'^,x'^) = {-{x3+X4,),-X2,-{xi+X4),Xi). 


( 6 . 2 ) 









FIG. 19: The xyyx interference contribution of figs. 5 and 6 with individual splitting times labeled. 
The numbers depict our labeling of the 4-particle Xj in (6.1). 


In this notation, one of the changes in going from xyyx to xyyx can then be summarized as 

'T* . \ 

•Jy 

Another difference between xyyx and xyyx can be seen by considering the splittings that 
occur at the start (fy) and end {ty) of the 4-particle evolution. At t=ty in £g. 19, particles 3 
and 4 are spectators to the splitting, and so the natural variables for describing the system 
(following our discussion in section IV B) are {C^^, 0 ^ 2 )- Similarly, the natural variables at 
the end of the 4-particle evolution, t=ty, are O^s). Looking only at the variable names, 
this seems like the same choices as for the xyyx analysis in (4.10). The difference is that 
(C'f4, Cl,) are the variables at the start of the 4-particle evolution of xyyx but the end of 
the 4-particle evolution of xyyx. As a result, Cl,) will be tied by the splitting matrix 
element to the initial 3-particle evolution {t^ < t < ty) for xyyx but was tied to the final 
3-particle evolution {ty < t < t^) for xyyx. The treatment of the initial and final 3-particle 
evolution was fairly symmetric in section V, and the effect of the 3-particle evolutions in 
the final result (5.45) appeared only through the parameters {Mi, 0.1) and (Mf,r2f). So one 
might expect that one can account for this difference between xyyx and xyyx simply by 
interchanging those parameters: {Mi,fli) y-)- {M{,Oi). In fact, this change is automatic 
under Xi — >■ x{ if one chooses to write 3-particle M and G in terms of the 4-particle Xi, as in 
the middle expressions in (5.2) and (5.4): 

Ml = xiXi{xi+X4)E —> x'ix'^{x'i+x'^)E = {xz+Xi)xiXzE = Mf, (6.3) 

Xi^x) 

and similarly Mf —)■ Mi and Gj -f-)- Gf. 

However, there is one additional element to the relation between xyyx and xyyx, which 
are the contractions of the helicity amplitudes associated with each splitting matrix to make 
the combination (4.38), 

a{x, y) + (3{x, y) + 7 (x, y) (6.4) 

discussed for xyyx back in section IV F, where the indices {m,n,m,n) were associated with 
the splitting vertices as in fig. 20a. Holding x and y fixed in a, (3, and 7 , the rules that 
we have discussed so far for relating xyyx to xyyx, which swap the initial and final 3- 
particle evolutions, would correspond to (6.4) contracted with vertices for xyyx as in fig. 
20b. However, as discussed earlier and exemplified by (4.37), the helicity sums should be 
the same for all of the diagrams in the hrst line of hgs. 5 and 6 . Reading from (4.37), the 
vertices should be labeled as in £g. 20c instead of £g. 20b. That is, we need to swap m n. 
From (6.4), this is equivalent to instead swapping a (3. 
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FIG. 20; Relation between labeling of splitting vertices and (6.4) for (a) our previous xyyx calcu¬ 
lation, (b) what we would get attempting to relate that calculation to xyyx by only taking Xj ^ 
[and so also (M, fl)i -H- (M, n)f], and (c) how the labeling should actually be for xyyx. 


Readers who find the above arguments a little too impressionistic should refer to appendix 
E. In particular, the discussion above completely sweeps under the rug a subtlety concerning 
the overall sign of the xyyx contribution, which is related to why we presciently replaced 
1 x 1 X 2 X 3 X 41 by —X 1 X 2 X 3 X 4 back in (5.39). (Though these are the same for xyyx, they differ 
by a sign for xyyx.) 

In summary, the steps necessary to change the result (5.45) for xyyx to a result for xyyx 
are 


• Xj —>■ x' (to include the formulas for Oy, and Oy); 

• (Mi, Gi) y-)- (Mf, Gf) if this change was overlooked when applying the Xi —)■ x' rule; 

• a{x,y) yy I3{x,y) . 


The interpretation of the 4-particle evolution time At in the resulting formula is now At = 
ty — ty rather than At = ty — ty, so that At is still positive. 

Applying these rules on the 1 /At divergence of (5.46), the corresponding divergence for 
xyyx is the At integral of 


iCl 




(Gf sgn Mf -|- Gi sgn Mj) x'ix' 2 x'^x'^ (Xf -f x^)^ (xg -t- X 4 ) 


/ 2 , 


At 


(/3 -\r Oi) 


Using (6.2), this can be rewritten as 


(/3 + 7)x^2x; 

(x'i-l-X4)(x3-fX4) 

(6.5) 


_(„ + ;3) _ (£±A2 £i', 

X1X3 

( 6 . 6 ) 

Note that the {a+f3) term will cancel when (5.46) and ( 6 . 6 ) are added together. To see the 
complete cancellation of all 1 /At divergences, we must go on to analyze other interference 
contributions. 


iCp^ai 
Ibvr^ At 


(Gi sgn Mi -|- Gf sgn Mf)x4X2X3X4(xi + X4)^(x3 + X4)" 


B. xyxy interference 

The xyxy interference is our name for the third diagram in figs. 5 and 6 . We will label 
the particles during the 4-particle evolution as in fig. 21 and designate the 4-particle Xj with 
tildes: 

(xi,X2,X3,X4) = (-?/,-(l- 2 /),x, 1-x). (6.7) 
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FIG. 21: The xyxy interference contribution of figs. 5 and 6 with individual splitting times labeled. 
The numbers depict our labeling of the 4-particle Xi as in (6.7). 


This is just a permutation of the x[ in (6.1), 


(xi, xs, xs, X4) = (x'a, x'^). (6.8) 

It is a convenient permutation of the labels because, as seen from comparing hgs. 19 and 
21 , it preserves the labels of which particles are involved in the splitting at the beginning 
and end of the 4-particle evolution: particles 3 and 4 are the spectators at the beginning 
(t=ty), and 1 and 3 are the spectators at the end {t=ty for xyyx and t=t^ for xyxy). So 
the structure of the two diagrams is the same, and the hrst rule of converting the results for 
xyyx is simply x[ Xi, or equivalently 1 O 2 and 3 O 4. Note that if we also apply this 
rule to the formulas for Mi and Mf in terms of x[, it gives the appropriate initial and hnal 
3-particle masses M associated with fig. 21: M, = {x'^+x'^)x' 4 ^x'^E does not change under 
13 ■(-)■ 24, and 

Mf = x[x'4^{x[+x'4)E -> Mf = XiX4{xi+X4)E = X2 X^^(x2Ex'^)E. (6.9) 

Here, 

Mf = —y{l—x){l—x—y)E (6.10) 

is different from our previous M’s because the particles in the final 3-particl_e evolution of 
xyxy are different than those for xyyx and xyyx. Unlike Mi and Mf, this Mf is negative. 
Similarly, the transformation gives us the appropriate 

Gf = y/ 

Unlike the Ui and Uf discussed previously, Gf is proportional to \/+z rather than \/^ and 
so has a positive rather than negative imaginary part. This is the reason that we earlier 
made the (at that time unnecessary) generalization from (5.8) to (5.9) in our discussion of 
integrating out the hrst and last splitting times. 

We can now express the transformation directly from xyyx to xyxy in one step by rewrit¬ 
ing ( 6 . 8 ) in terms of the xp 


% n ^ 1_ l_\ ^ I tqA f I ^ I I \ 

2E yxf X4 X1+X4J y 2E \ y 1—x 1—x—yJ 


(Xf, X2, X3, X4) = (-X2, -(X3+X4), X4, -(xi-t-X4)) . (6.12) 

Now remember that, by itself, the transformation from x* to x' gave us hg. 20b instead of 
what we needed for xyyx, hg. 20c. If we now make the further transformation 12 y-)- 34 
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FIG. 22: As figs. 20(b,c) but here for xyxy. 


that takes us from x' to Xj, fig. 20b will become fig. 22b for xyxy.^^ But what we need, as 
read from (4.37), is £g. 22c. This requires mnn —)■ nnm, which, from (6.4), is equivalent to 
af3'y qtt/d. 

In summary, the steps necessary to change the result (5.45) for xyyx to a result for xyxy 
are 

• Xi ^ Xi (to include the formulas for r 2 ±, Oy, and Oy); 

• (Mi, Gi; Mf, Gf) (Mf, Gf; Mi, Gi) if this change was overlooked when applying the 
Xi —)■ Xj rule; 

• {a{x,y),l3{x,y),'^{x,y)) ( 7 (x, ?/), a(x, ?/),/3(x, ?/)) . 

The interpretation of the 4-particle evolution time At in the resulting formula is now Af = 

0 . 

Applying these rules on the 1/At divergence of (5.46), the corresponding divergence for 
xyyx is the At integral of 



(7 + /?)x2X4 


(Xi-f-X4)(x3-fX4) 

(6.13) 


Using (6.12), this can be rewritten as 



(X 1 +X 4 )(X 3 - 1 -X 4 ) X 1 X 3 


Note that all of the Gi terms now cancel between (5.46), ( 6 . 6 ), and (6.14)—that is, between 
the three diagrams shown in the hrst line of hgs. 5 and 6 . 

VII. 1/At DIVERGENCES AND ie PRESCRIPTIONS 

We have now seen a subset of 1/At divergences cancel in the hrst line of kg. 6 . To get 
the rest of the cancellations just requires a slightly bigger subset of diagrams: All of the 


Again, see Appendix E is this be less than clear. 
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yxxy yxxy yxyx 

FIG. 23: A subset of crossed diagrams for which all 1/At divergences cancel. 


1/At divergences cancel between the six diagrams shown in £g. 23. The second line of fig. 
23 is just the conjugate of the hrst line permuted hy x ^ y. One may check that adding 
together the 1/At terms of (5.46), (6.6), (6.14) and their conjugates with x -H- y gives zero. 

The complete set of crossed diagrams depicted by £g. 6 can be written as the sum of £g. 
23 and its conjugate, plus the remaining possible permutations of the three hnal daughters 
{x,y, 1—x—y). So 1/At divergences will cancel for the complete set of crossed diagrams as 
well. That means that, if we combine results for the crossed diagrams before doing the At 
integral, and we subtract the corresponding vacuum results as in (5.47) and (5.48), the At 
integral will be hnite and can be performed numerically. Before we summarize formulas for 
the hnal result, however, we need to take care of a subtlety with the contribution from At=0 
to which we have previously alluded. And before we do that, it will be enlightening to first 
review some technical details of the calculation of single (as opposed to double) splitting. 

Note added'. The analysis of this section misses some terms in the contribution from 
At = 0. A corrected analysis, which includes the missing terms, may be found in ref. [32]. 


A. ie prescription in single splitting 


Return to the formula (2.37) for the single splitting rate. Using the propagator (2.38), 
this formula becomes 

^ ^ _ aPi^3{x) r ^2 At), (7.1) 

dx 71 Jq 

where At = t^ — t^- The integrand blows up like l/(Af)^ as At —)• 0. A simple way to 
avoid this problem altogether is to subtract out the vanishing rate of splitting in vacuum, 
analogous to (5.47): 


dr _ dr 

dx dx 


[dr 


vacuum 


dx 


dv 

dx dx 



rig csc^(r2o At) 


1 

AT. 


(7.2) 


This integral is convergent and gives the result 


dr 

dx 


aPi^3{x) 

71 


Re(tf2o), 


(7.3) 
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which in turn gives (1.5a). 

However, it will be useful to see how to get the result without making the vacuum 
subtraction. We’ll see how to deal with the l/(At)^ singularity in (7.1) directly, but the 
argument will perhaps look more familiar if we hrst transform the At integral from 0 to cx) 
to an integral from — cx) to +cxd. 


1. An integral from —oo to +oo 

For the sake of interpretation of At, it will be useful to consider how we would have 
arrived at (7.1) if we had computed the xx interference diagram of £g. 10 and its conjugate 
XX separately, as depicted in £g. 24. The 3-particle given by (2.33b), 


n 


2E \xi X 2 X 3 J 


(7.4) 


gives (1.5b), 



( 7 . 5 ) 


for the XX interference diagram [for which {xi,X 2 , X3) = (—1, 1—x, x)]. But it gives fig for the 
XX diagram [for which [xi,X2,X3) = (l, —(1—x), —x)]. The results for these two diagrams, 
after following the same steps as in section II, are 


'dV 

dx 


XX 


aPi^3{x) 

2ti 


d{At) fig csc^(f2o At) 


with At = t^ — fx and 


'dV 

dx 


XX 


aP] 


l^>3 


X 


2tt 


d{At) fI*oCsc2(fIo At) 


(7.6a) 


(7.6b) 


with At = p — tx- Adding them together reproduces (7.1), as it should. Now let’s take 
At -k- —At in (7.6b) to write the sum of (7.6) as 


dr _ aPi^sjx) At), At<0-, 

dx 271 l^flg CSC^(ng At), At > 0, 


where now At = tx — tx in both cases. 

If the integrand in (7.7) were an analytic function of At, we could now say that the 
problem of the singularity at At = 0 is just a problem of understanding how one should 
integrate around the pole. The integrand is not analytic because of the various high-energy 
approximations that have been made, but the divergent piece at At = 0 is analytic. We can 
rewrite (7.7) as the sum of a convergent integral 


^ ^ aPi^sjx) /•+“ IQf csc^Q*o At) - (At)-^, At < 0; 
dx 271 csc^(fto At) — (At)“^, At > 0 
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FIG. 24: The interference diagrams contributing to single splitting. 


[which is equivalent to (7.2)] and a potentially divergent integral 

'dri /■+“ d{At) 


(7.9) 


It doesn’t really matter how we integrate around the double pole in (7.9) because we will 
get zero however we do it. However, getting straight the prescription will help deal with the 
1/At divergences that arise in the double splitting case. So we now discuss the ie prescription 
that should be used for At. 


2. The ie prescription 


Consider the 6H matrix elements that appeared in the formula (2.2) for single splitting, 


{Pfki,t^\p^ki,Q{p^k{\-i6H\p/) \{ / {pi,t^\p^,t^){p^\i6H\pfkf) ]. (7.10) 

'pi,ki \Jpx J \Jpx 


Using completeness relations, (7.10) can be rewritten as 


(7.11) 


This is a Wightman correlator: the SH for the splitting at t^ is always ordered to the left of 
the SH for the splitting at regardless of the time order of H and H. The ie prescription 
for the small-At divergence of a Wightman correlator is 


t—t' —)■ t—t'—ie for {A{t) B{t')). (7-12) 

See, for example, the argument in Wilson [25] with regard to the singularity behavior of 
A{x) B{y) in the operator product expansion. (We provide our own summary of the argu¬ 
ment in appendix A.) In our context, we can summarize this rule as follows: 

Regard times in the conjugate amplitude as displaced by —ie with respect to 
times in the amplitude. 


3. Other ways to use the prescription 

In our analysis of double splitting, we have always dehned our At to be positive. In 
single splitting, the analog is to use the 0 < At < oo integrals of (7.1) or (7.6) instead of the 
—oo < At < oo integral of (7.7). It is useful to notice that the ie prescription works just as 
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well in this form. In the case of (7.1), At = — t^. Since is in the amplitude and in 

the conjugate amplitude, (7.1) becomes 


dr aPi^3{x) 


dx 71 

The divergent piece is proportional to 

poo 

Re 


Re / d(At) hip csc^(r2o (A^ — *e)). 


( 7 . 13 ) 


=Re^=0, 


'll - «) 


-Z€ 


(7.14) 


Equivalently (and most like what we will do for double splitting) we can get the same 
result from (7.6) by realizing that At = ^ At — ie in the xx piece (7.6a) but 

At = tx — t At + ie in the xx piece (7.6b), so that the divergent piece of the sum is 
proportional to 

r d(At) , r 

I (At+tef ' > 


B. Consequences for double splitting 

Just like in the single splitting case, the vacuum l/(At)^ divergences in double splitting 
will be canceled by adding diagrams to their conjugates, or can be discarded by simply 
subtracting the (vanishing) vacuum result. However, the sub-leading 1/At divergences give 
something non-trivial. 


1. A warm-up example 


As a simple, warm-up example from calculus, imagine the integral 


I = 



n 

At — ie 


iTrhl, 


( 7 . 16 ) 


for some complex frequency H and some bound R. We may replace At ^ —At for At < 0 
to rewrite the integral as 


I 



At — ie~^ 



-n 

At ie’’ 


(7.17) 


which is a 1/At analogy to the l/(At)^ integrals we had in (7.15). Now combine the two 
integrands in (7.17). If we ignored the important ie prescriptions in (7.17), we would naively 
conclude that the 1/At terms cancel each other, leading to the incorrect answer 1 = 0. 

One can equivalently write (7.17) as 


I 



4(A1)A + 



d(At) 


-O 


( 7 . 18 ) 


where the contours Ci and C 2 are shown in Eg. 25. The integrals indeed cancel except for the 
contributions from the tiny region around At = 0, where the contours differ. The integrals 
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FIG. 25; The contours Ci and C 2 of (7.18) in the complex At plane. 


around each quarter-circle give ±i 7 r /2 times the residue at the origin, which brings (7.17) 
to the same result I = inQ as (7.16). 

A quick mnemonic for all of this is to use the standard relation 


1 

At =F 



± in 6{At), 


(7.19a) 


where “P. P.” indicates the principal part prescription, and with the understanding that 

R 

d{At)5{At) = \ (7.19b) 

integrates over only half of the 5 function. The principal part pieces then cancel between 
the two terms in (7.17), and the 5 function pieces account for the answer I = inVL. (In 
this context, the principal part prescription does not really do anything since we are not 
integrating over negative At.) 



2. Double splitting 


Let us now apply (7.19) to the 1/At divergences found previously for double splitting. 
For the xyyx result, this divergence was given in (5.46), where At = ty — ty. Since times 
in the conjugate amplitude should get a —ie, we take At ^ At — ie in (5.46). After then 
applying (7.19), we will have a 1/At piece and a +in6{At) piece. The 1/At pieces will 
cancel between diagrams, as already noted in the discussion of £g. 23. The +in6{At) piece 
will give an additional contribution to the xyyx interference that we need to retain, given 
by 


dT 
dx dy 


pole 


xyyx 


cWs 

32n 


(Oi sgn Mi -|- Of sgn Mf^xlx 2 xlx 4 {xi + Xi)‘^{x^ £ 4 )^ 


X 


{a. fd) 


{a + 7)52^4 

{XI+X4){X3+X4) 


(7.20) 


For xyyx, we have a similar situation, but At in ( 6 . 6 ) is At = ty — ty —)■ At -|- ie, and 
so the inS{At) in (7.19) has the opposite sign. As a result, the corresponding contribution 
arising from ( 6 . 6 ) is 


dT 

dx dy 


pole 


xyyx 


cWs 

32n 


(Oi sgn Mi -|- Of sgn Mf^xlx2xlx4{xi + Xi)‘^{x 3 -1- x^Y 


X 


{a. -|- /3) 


{(3 + '^)x2Xa 


X1X3 


(7.21) 
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FIG. 26: As fig. 18 but instead for the xyxy interference of fig. 21. 


The xyxy contribution of fig. 21 is a little more subtle because there At = — ty does 

not pick up a Aie prescription from what we have discussed so far. Here, we will indicate 
how one might guess the answer, and we leave a more precise argument to appendix D3. 
Analogous to the earlier discussion of fig. 18 for xyyx, the VtJAt and Hf/At divergences 
of (6.14) arise from the situations depicted in fig. 26. Looking at the figure for the Hf/At 
divergence, remember that, in the derivations of our expressions (section VA in particular), 
we have already integrated the first time t^ all the way up to the second time (in this case 
ty). And so one might guess that the t^ — ty will inherit the ie prescription of t^ — t^, which 
is tx — tx — Our guess is that At = t^ — ty ^ At — ie will give the correct result for the 
Clf/At divergence. 

In contrast, for the VL-JAt divergence in fig. 26, t^ is not close to ty, but instead ty is 
close to tx. One might guess that t^ — ty then inherits the ie prescription of ty — ty, which 
is ty — ty + ie. So our guess is that At ^ At + ie will give the correct result for the fli/At 
divergence of xyxy. These guesses are correct (see appendix D 3) and put together with 
(6.14) and (7.19) give 


dr 

dx dy 


pole (J2 ^,2 


A Oi ~ ~ 

(OiSgnMi - VLiSgnMi)xlx2xlxi{xi +Xif{x^ + Xi) 


xyxy 


X 


{a + 'y)x2XA + l)x2Xi 


_{X1+X4){X3+X4) 


X1X3 


(7.22) 


The sum of the 0; terms of (7.20-7.22) is 


IQtt 


ilixlx2X3X4,{xi + Xi)‘^{x3 + £ 4 )" 


— (a + /?) + 


{a + 'y)x2X4: 

{XI+X4){X3+X4) 


(7.23) 


(where we’ve finally used the fact that sgnMi = +1 since we will no longer be considering 
variable substitutions that will change the sign). If we add (7.20-7.22) to their conjugates 
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with X y, as in the subset of diagrams of £g. 23, then we get 


dT 

dx dy 


pole 
fig. 23 


CM 

IhTT 


(Oj + Of) a;fa; 2 a; 3 a; 4 (a;i + £ 4)^(£3 + £ 4 )" 


X 


— (q; + /3) + 


{a + 'y)x2X4: 

{X1+X4){X3+X4) 


+ {x y}*. (7.24) 


(where we have now used the fact that sguMf = — 1 ). So, even though the l/(At) terms 
naively cancel between the diagrams of fig. 23, the pole contributions arising from the ±ie’s 
do not. 

Recalling that the (q;,/ 5 , 7 ) of (4.39) are symmetric under x -H- y, (7.24) can be written 
more explicitly as 


~ dr 1 

_dx dy_ gg 23 


where 




A«s 


^-(i-x)^l-x-y,y\ 

a + ^^xy 


X xy{l-x) (1-1/) {l-x-yY 


(cr + /3) + 


(l-a;)(l-//) 


(7.25a) 


O 


tCi,X2,X3 



(- + - + -). 

\Xl X2 X3J 


(7.25b) 


VIII. SUMMARY OF RESULT 


We now collect our final result for the sum of all crossed diagrams, depicted by fig. 5 or 6 , 
for gluons in the large-iVc limit in a form appropriate for (convergent) numerics. The result 

is 


dr 


dx dy 


crossed 


A(x,y) + A(l-x-p,p) + A(x, l-x-p), 


( 8 . 1 ) 


where A(x, y) is the result for the six diagrams of fig. 23 plus their conjugates, with vacuum 
result subtracted. We will write this as 


A{x,y) = 2Re 
= 2Re 


dV 

dx dy 


fig. 23 


lim 


dV 

dx dy 


J fig. 23 


dr 

dx dy 


pole 


r*H-oo 


+ 


J fig. 23 


d(At) 2 Re \B {x, y, At) + B(y, x, At)], ( 8 . 2 ) 


where B [x, y) is the At integrand for the first three diagrams of fig. 23 (with vacuum results 
subtracted). This is 

B{x, y, At) = C{{xi}, a, fd, 7 , + C{{x'i}, fd, a, 7 , At) C{{xi}, 7 , a, (d, At) 

= C(-l, y, l-x-y, X, a. Id, 7 , At) C{-{l-y), -y, 1-x, x, Id, a, 7 , At) 

+ Ci-y,-{^-y),x, l-x, 7 , a. Id, At), (8.3) 
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where the three terms correspond to xyyx, xyyx, and xyxy respectively, and a{x, y), 13{x, y), 
and 7 ( 0 ;, y) are given by (4.39). C is dehned to have the vacuum result subtracted, so write 

C = D-\imD. (8.4) 

D, dehned as the integrand for xyyx, is given by (5.45) as 


D{xi,X 2, X3, X4,, a, ( 3 , 7 , At) = 

^ {-xiX 2 X 3 XA)^+Vt_ csc(n+At) csc(n_ At) 

^ + aYyyYyy)I„ Y {tt Y ^ Y 20')Zyy/l 

+ [{aYY)YyY,Y{DYY)Yy,Yyy]h-(aYf>YY)(Yy,YyhYYyYyyU)] (8.5) 

with the In dehned by (5.44); the (A, Y, Z)’s dehned by (5.42); the 4-particle normal modes 
and frequencies used in their dehnition given by (5.21-5.25), (5.28), (5.33), and (5.40); and 
the 3-particle M’s and ft’s in these formulas dehned specihcally in terms of the 4-particle Xi 
(the arguments of D above) as 


Mi = xiXi{xi+X4)E, 


Mf = X3XA{x3+Xi)E 


and 


( 8 . 6 ) 





— + —- —) 

Xi X 4 X1+X4 J 





—- 

\X3 X4 X3+X4J 


(8.7) 


The pole contribution in (8.2) is given by (7.25). 

The g —)■ 0 limit for the vacuum piece in (8.4) corresponds to taking all fTs to zero. The 
vacuum piece in (8.4) could be worked out algebraically, but it is simpler to just use the 
same numerical code as for (8.5) but with the replacements 


hli —)■ 0, hlf —)■ 0, ficot(fi At) —)■ (At) ficsc(fi At) (At) ^ ( 8 . 8 ) 


[and so also r2± csc(r2±At) (At)“^ in the prefactor of (8.5)]. 

It is possible to scale out the factors of qa and E from the above expressions by replacing 
At by the dimensionless variable At = {qa/EY^'^A t. For numerics, it is convenient to work 
in units where qA=3 and E=l, which then gives the result for the rate dV/dxdy in units of 
(to/S)'-'". 


IX. BEHAVIOR OF RESULT 

Though this paper calculates only a subset of the diagrams contributing to dV/dxdy, it 
is interesting to examine the behavior of the results in the soft bremsstrahlung limit. 
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A. Comparison to earlier results for y x \ 

As mentioned before, refs. [ 6 - 8 ] have studied the effect of double splitting on energy loss 
when one of two soft bremsstrahlung gluons becomes extremely soft compared to the other 
(y a; 1). They hnd that the double-splitting correction to energy loss is enhanced 
by a double logarithm that, at leading logarithmic order, can be absorbed into the single 
splitting result (1.5) for energy loss by a redehnition of q. This energy-dependent correction 
5q to q is seemingly universal and was computed yet earlier by Liou, Mueller, and Wu [11], 
with result 

(9-1) 

In the current context of a thick medium, the relevant distance scale L is the formation time 
associated with emission of the x gluon; 



[The fact that this is a parametric statement (~) rather than an exact equality does not 
matter at the level of identifying the leading logarithmic result.] The /q above characterizes 
the small distance scale, characteristic of the medium, at which the multiple scattering 
(harmonic oscillator) approximation breaks down. A discussion of resumming the leading 
logarithmic terms at all orders may be found in refs. [ 6 , 8 , 11, 12]. The results are very 
interesting because they can change the power law dependence of how energy loss depends 
on energy. For instance, the stopping distance of a high-energy parton due to repeated 
single-splitting processes scales with energy as but Blaizot and Mehtar-Tani [ 6 ] have 
discussed how resumming the double logarithm above can modify the scaling to 

These earlier results provide a useful check on our own results: Do they agree if we 
take the y a; 1 limit? We will focus on a comparison with Blaizot and Mehtar-Tani 
[ 6 ] and Wu [ 8 ]. Comparison is slightly complicated by the fact that in this paper we have 
only calculated certain contributions to dT/dxdy but have not calculated virtual corrections 
like fig. 8 that are needed for the energy loss calculation. On the flip side, the processes 
studied by refs. [ 6 , 8 ] do not include all possible time orderings of the splittings (presumably 
because only certain orderings are relevant to the double logarithm). So, we will need to 
compare a subset of our diagrams to a subset of theirs. In our terminology, this subset is 
xyyx + xyyx + zyyz + zyyz (plus their conjugates). For reasons explained in Appendix F3, 
the xyyx -\- xyyx contributions are difficult to disentangle from the virtual processes also 
included in the analysis of refs. [ 6 , 8 ], and so we shall content ourselves with checking the 
contribution from zyyz + zyyz, which are just the xyyx and xyyx diagrams analyzed earlier 
with the substitution x ^ z = 1—x—y. An alternative depiction of zyyz was given earlier 
in £g. 9. (In Wu [ 8 ], our zyyz and zyyz are respectively the second diagram of his (30) and 
the first diagram of his (32).^^ In Blaizot and Mehtar-Tani [ 6 ], they are called A 3 and Ci.) 

The work of Wu [ 8 ] is closest in technique to our own. In appendix F, we briefly sketch the 
relationship between Wu’s expressions and the ?/ < a; 1 limit of our formal expressions. 


Here and throughout, we refer to Wu’s equations [8] by their numbering in the arXiv version of the paper 
(versions 1 or 2) rather than the JHEP version of the paper, in part to make it a little easier to distinguish 
references to our equations from references to his. 
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Wu then goes on to extract the double log behavior (getting the same result as others) to 
which we now turn. 

The contributions from zyyz + zyyz that give rise to double log corrections turn out to 


be 



^«A*« 



(9.3) 


which is equivalent to the regions found to give the double logarithms in refs. [6, 8] (see 
appendix F for details). When we expand our own final results (from section VIII) in this 
range, also taking y a; 1, we find 


2 Re 


dV 


dx dy 

zyyz+zyyz 


(2 + 0 + 0 ) 


1 j (12-3 + 1) 


n‘^xy{Aty 


+ 


Sir'^yx^/^At 


+ ... 


Ciat 



(9.4) 

Here, we have distinguished which contributions arise from interaction of the y gluon with 
the medium, using notation that we will explain in a moment. Specihcally, the evolution of 
the y gluon for zyyz + zyyz is only relevant during the 4-particle part of the evolution, during 
which non-zero q only enters the calculation via non-zero values of In the y a; 1 
limit, the relevant values of (5.21) approach 






-^qA 

2yE' 




-3igA 

8xE 


(9.5) 


for zyyz and 

—y —y —(9-6) 

for zyyz. The “(12 — 3 + 1)” in (9.4) means that the contributions that do not depend at all 
on hlj- (and so correspond to no interactions of the y gauge boson with the medium) give 
“12” there; contributions that are proportional to give —3; and contributions that are 
proportional to give +1. Similarly for the “(2+0+0).” 

The l/(At)^ and l/(At) pieces of (9.4) have nothing to do with the small-At divergences 
discussed back in sections VD and VII. Those divergences and their cancellations arise 
in the limit of arbitrarily small At rather than the At range (9.3) being considered here. 
(Specifically, the small-At cancellations in [zyyz + zyyz + zyzy] + [z y]* require At 
^/yJx\/yE/q.) 

The double log effects which are found in the energy loss calculations of refs. [6, 8] arise 
from single interactions of the soft y gluons with the medium, and so, when comparing 
similar diagrams, correspond to just the pieces 


(0 3 + 1) 2 2 

M^yx^At^^ E 


(9.7) 


of (9.4). Now integrate over y and At, and compare to the single splitting result (1.5) for 
dr/dx in the limit x 1, which is 


dr CAtts [Ia 

_ _ A / _ 

dx TTo;^/^ V E 
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The result is that the double-splitting contribution considered above can be absorbed (up 
to higher order corrections) into the single splitting result (9.8) by ija —t (?a + Sqa with (see 
Appendix A) 


e.. (~3-|-l) GsCa ^ 

SqA ^ ^- qa 


TT 


dyd{At) (-3-hl) OsCa . 2 f 


(9.9) 


[where now we show only the contributions arising from and VLy in the notation 
“(—3-1-1)”]. Even after adding in xyyx + xyyx (and conjugate) contributions, this will 
not exactly match (9.1) because in this paper we have not included the virtual corrections 
relevant to energy loss. However, we have checked that the contribution to Sqa from just 
the zyyz + zyyz (and conjugate) diagrams in refs. [6, 8] would have led to the same result 
for the double log pieces, i.e. 


dqA _ (-3+1) 5qA 

qA \ ^yyz+^yyz 4 q^ 

+conjugates 


(9.10) 


with dqAjqA on the right-hand side representing the total 6qA/qA — («s/27r) ln^(L//o) of 
(9.1). 


B. Small y behavior of full result 


Our calculation contains other diagrams and effects besides the ones considered in the 
double log works of [6, 8]. Numerically taking the y a; 1 limit of our full result for the 
crossed diagram contributions to dT/dxdy, we hnd that the result behaves as^^ 


dT 

dx dy 


oc 


crossed 


Mv/x) 

xy^A ’ 


(9.11) 


as shown, for example, by the numerical results of £g. 27.^^ The crossed contribution to 
dr /dx dy is mildly divergent in the infrared. If we cut off x and y at some lower value 
a^min, (9-11) implies that the crossed contribution to the total rate T for double splitting 
into non-virtual daughters scales as In Xmm- For comparison, the total rate for single 
splitting diverges as a relatively mild x ^J^ . 

Keep in mind that we have only computed a subset of diagrams (fig. 6). Our crossed 
interference diagrams do not sum to be the square of something and so, unlike total dT/dx dy, 
the result for [dT/da; dy]crossed need not be positive. Also, with regard to applications to 
energy loss, remember that we have not included virtual corrections to single bremsstrahlung. 


See Appendix A for a discussion of how the pole terms (7.25) make a partial but incomplete cancellation 
of the leading behavior of (9.11). 

As an additional numerical check for anyone who ever desires to implement our formulas of section Vlll, 
our value of [dr/dx dy]crossed is 11.463126196457 C\q^\/ q^lE for {x,y) = (0.3, 0.6). [Note added. The 
correction to the pole terms found in ref. [32] changes this to —10.892657927744 C’^a'^y^qA/E-] 
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FIG. 27: Our numerical results (solid lines) for the total crossed diagram contribution to dV/dxdy 
[in units of qx/E] vs. y for fixed x = 10“^. The solid violet and blue curves depict 

positive and negative numerical values, respectively. The dT oc hi{y/x) dashed line shows the 
y-3/2 in(y/a:) behavior of the y <C a: <C 1 power law quoted in (9.11). The dT oc y~^ ln(x/y) dashed 
line shows the x~^ ln(x/y) behavior of the same power law if one switches the labels x and y. We 
have only shown results for y < 0.5; results for y > 0.5 are given by the permutation symmetry 
y ^ z = 1—x—y of the problem. [Note added. This plot does not reflect the correction to the pole 
terms found in ref. [32], which change the result to be everywhere negative.] 


X. CONCLUSION 

This paper is long. Congratulations if you read (or even skimmed) this far! 

A complete calculation of dT /dx dy outside of the x -C 1 and/or double log approximation 
will require a little further work. But the analysis in this paper of all of the crossed diagrams 
provides most of the techniques that will be necessary for the full analysis, and we look 
forward to filling in the rest. 

The reason that we have saved un-crossed diagrams such as fig. 7 for later work is that 
their 1/At divergences do not cancel as simply as fig. 23. One must instead take great care 
in separating the overlapping formation time corrections for double bremsstrahlung from 
the naive calculation of two independent, consecutive single-bremsstrahlung processes. But 
that is a story for another day. 


Acknowledgment s 

This work was supported, in part, by the U.S. Department of Energy under Grant No. DE- 
SC0007984. We thank Han-Chih Chang, Gabriel Wong, and Bin Wu for useful discussions. 


61 











Appendix A: More details on some formnlas 


Eq. (2.1) and (2.2): We use non-relativistic normalization convention for \p), so that 
{p'\p) = — p') and f = f d^p/{2nY. The factor of is just the square 

of the normalization factor associated with the initial state |pi), since {Pi\Pi) = V± in our 
normalization. [In these formulas, we have been a little sloppy about factors of V± associated 
with the normalization of the photon state because we wanted to emphasize that fcf = k^. 
If concerned about the details of all factors of V±, refer to the more general formula (2.34) 
instead. If you want to see how engineering dimensions match up, then also see (4.29-4.31) 
and the discussion of them in this appendix.] The factor of E/2'k in (2.1) comes from 
rewriting the longitudinal phase space factor dkz/27i in 


dr 


d^k 

(^ ■■■ 


(Al) 


as dkz/271 = {E/27 i) dx and dividing both sides by dx. Evolution of the hnal state \pikf) 
after can be ignored by unitarity of time evolution (before taking the statistical average 
((■ ■ ■))) and the completeness relation |pffcf)(pffcf| = 1 within the two-particle sector. 
See section IV A for a more detailed argument in a different but closely related context. (The 
same hidden assumption discussed in footnote 21 applies to the argument here.) A similar 
argument can be made for time evolution of the initial state before by invoking rotational 
invariance. Rotational invariance of the rate we are computing means that we can average 
the initial Pi over small rotations, which in the context of our large-pz approximation is 
discussed in detail in appendix B 2, with the result that one may replace |pi)(pi| by something 
proportional to |Pi)(Pi|- The same unitarity and completeness arguments used for the 
hnal state for t > then imply that we also need not follow the evolution of the initial state 
for t < in ( 2 . 1 ) and ( 2 . 2 ). 

Eq. (2.21): Here is a more general argument for (2.21). In general, any translation- 
invariant harmonic approximation for 3 particles can be put into the form 


^ (^ 1 ) ^ 2 ; ha) — A 2 i(h 2 ~ hi)^ -|- ^ 32(^3 — h 2 )^ + ^ 13(^1 — (A2) 


with some coefficients A^j. In the special case of bi = b 2 , this gives 


y (hi, hi, ha) — (A32 -|- Ai3)(h3 — hi)^. 


(A3) 


However, in this case we have color charge T 3 at ha and total color charge Ti -|- T 2 = —T 3 
at hi = h 2 - This is effectively a two-particle problem with separation b = b^ — bi, and the 
small-h limit of that problem is precisely what is used in proposals to non-perturbatively 
dehne (using real-time Wilson loops [21]) the parameter q for color representation T 3 , which 
means that 

V(hi,hi,h3) = -ig3(h3-hi)2. (A4) 

Combined with (A3), this gives a constraint on the values of the coefficients Aij. Permuting 
the particle labels on this argument gives three constraints on the three unknown coefficients 
Aij, which determines (A 2 ) to be ( 2 . 21 ). 
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Eq. (2.22): In weak coupling, q defined by (2.22) is UV logarithmically divergent, 
oc a^n f (Pq±/q]_ where n is the density of partons in the medium. In the harmonic approxi¬ 
mation, q then has a mild logarithmic dependence on h, and the harmonic approximation is 
effectively a leading-log approximation (to which sub-leading-log corrections may be system¬ 
atically computed^^). However, if the coupling is replaced by running as{q±), q defined 
by (2.22) is convergent. This form of q is only relevant if the total momentum transfer Q± 
during a formation time is large enough that the integral ( 2 . 22 ) is not significantly modihed 
if restricted to q± < Q±. This happens when q;s(Q±) cts^mj)), where mo is the inverse 
screening distance in the medium. For thick media, Q± ~ {qEY^^ for hard bremsstrahlung, 
and note that our assumption throughout this paper is that as(Q±) can be treated as small. 
Finally, qL describes the average Q\ picked up from the medium over a distance L and 
characterizes a Gaussian probability distribution for Q± after a large number of collisions. 
However, power-law tails of the probability distribution for Q±, representing rare collisions, 
can dominate certain types of calculations.^^ 

Eq. (4. 16): In this formula, we make use of = 0 to occasionally rewrite, for 

example, —£ 3—£4 as £i+£ 2 . Other than that, the expression follows from the rules of 
£g. 16 together with the overall factor of {E/2t:Y of (4.10), which comes from rewriting 
dkz/271 = {E/27 i) dx and dKzl 2 'n = {E/ 2 ti) dy, similar to the single factor of E/ 2 ti in ( 2 . 1 ). 
The one thing to be careful about in applying the rules is the signs of the Bij. For instance, 
consider the rule associated with the vertex at t = in fig. 15. Using the top rule in fig. 16, 
the vertex is associated with a factor of V5^^^(Bjj) according to the labeling of that rule, 
but we may just as well write V(5^^^(Bifc) using the identity Bji = Bkj = Bik noted in the 
caption. Comparing fig. 16 to fig. 15, this V(5*^^^(Bjfc) translates to V( 5 *^^^(.Bi 4 ) in the index 
labeling used in the latter figure. In (4.16), however, we have used C44 at the start t = ty of 
the 4-particle evolution. As a result of this choice, it was then convenient to define our B 
for the initial 3-particle evolution E < t < ty as B = B 41 instead of .B = B 14 . So our factor 
of V(5^^^(Bi4) at t = E is — V(5*^^^(B^). In contrast, the B defined for the final 3-particle 
evolution ty < t < is B = B 34 , which is the same B that comes from applying the top 
rule of fig. 16 to the splitting at t = t^. So the corresponding factor for that splitting is 
-1-V(5^^^(B^). In a similar vein, the rules of fig. 16 give factors of ^ 

at t = ty and at t = ty. We then integrate all the ^-functions (integrating by 

parts as necessary) to get (4.16). 

Eqs. ( 4 . 29 - 4 - 31 ): We offer here a quick note on the scaling dimensions of our formulas. 
Since we use non-relativistic normalization {p\p') = {27rY5^‘^Yp — p'), then |pi) has dimen¬ 
sions mass“^ and {Pj,Pk\ has dimension mass“^, and so one would expect {pj,Pk\SH\pi) to 
have dimension mass“^. The right-hand side of (4.29) has an implicit {2TiY5^‘^Ypi—pj —pk)- 
Including this factor gives the right-hand side dimensions of mass“®/^. The discrepancy 
has to do with the normalization of the states of the longitudinal momentum. Here 
are two ways to handle the normalization. Method 1: If we normalize states so that 
{p,x\p',x') = {27rY5^'^Yp — p')5x,x'-i with a Kronecker 5 for the x’s, then there is really 
a factor of 14 ^ on the right-hand side of (4.29), where I 4 is the (infinite) length of the 
; 2 -direction. (This I 4 should not be confused with the finite physical width of the medium 


See, for example, refs. [26, 27]. 

See, for example, the discussion of power law tails in sec. 3.1 of BDMPS [4]; Zakharov [28]; Arnold [29]; 
appendix A of Peigne and Smilga [30]; and D’Eramo et al. [31]. 
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in cases where the medium is hnite.) This fixes up the dimensions of (4.29). With this 
normalization convention for states of longitudinal momentum, the sum over the longitu¬ 
dinal momentum of the bremsstrahlung particle in single splitting, for example, becomes 
f Vzdpz/271 instead of f dpz/2n, and so the single-splitting formula ( 2 . 1 ) should have in- 
eluded an overall factor of 14- This 14 then cancels the two factors of 14 that we get 
from the two dH matrix elements. Similar cancellations occur in the double-splitting case. 
We have not bothered to include any of these canceling factors of 14 in our formulas. Method 
2 : Alternatively, one could normalize states as (p, a;|p', x') = (27r)^(5*^^^(p — p') 5{pz —p'^) and 
avoid ever introducing 14. This normalization would change the dimensional analysis to 

_1 /o 

make the two sides of (4.29) agree dimensionally, without any additional factor of 14 , if 

we take there to be an implicit {pi—pj—pk) d{piz—pjz—Pkz) on the right-hand side. 

Eq. (5.9): If hi has a negative imaginary part, we get cot(r2cxD) = i and (5.8). If hi 
had a positive imaginary part, we would instead get cot(r2cxD) = —i and so (5.8) with 
exp(—I replaced by exp(-|-| Mhli?^). We could therefore write the exponential in 

the general case as expQ sgn(Imhi)). However, there is a correlation between the 

sign of M and the sign of ImH that makes for a more compact expression. First note that 
if the 3-particle effective Hamiltonian EL of (2.33) makes any sense for our problem, then 
it should have a negative imaginary part, corresponding to exponential decay rather than 
exponential growth. This means it should have MImHQ < 0, and so, also using (2.33b), 

sgn Im H = sgn Im = — sgn M. (A5) 

This allows us to rewrite our general-case exponential above in the form of (5.9a), as 
promised. 

Readers may wonder whether the condition that MImHg < 0 could ever be violated. 
Using xi + X 2 + Q for 3-particle evolution, (2.27), and (2.33b), 

MlmH^ = -\{{q2 +q-i-qi)x{ + {q-i+qi-q2)xl + {qi+q2-q?)xfj. (A 6 ) 

Requiring the right-hand side to be negative for all allowed values of Xi satisfying xi+ X 2 + 
Xa = 0 is equivalent to requiring that the largest value of 1 /^ be smaller than the sum of 
the two others.^® So, if gs is the largest, the condition would be 

< \f4i + V^- (A7) 

In the weakly-coupled case, this would be < ^/Ci -|- which can be rewritten as 

IIT 1 +T 2 II < IIT 4 I + IIT 2 II, (A 8 ) 

where Ti are the color generators associated with each particle, acting on a 3-particle color 
singlet. The inequality (A 8 ) follows from the triangle inequality, given that the action of Ti 


One (perhaps inelegant) way to derive this condition is to plug = — (a;i -l-a; 2 ) into (A 6 ) and then rewrite 
the condition in terms of the ratio r = xi(x 2 , which gives q 2 r^ + ( 91 - 1 - 92 — 93 )^ -t 91 > 0 , which we require 
for all values of r. Requiring the minimum with respect to r to be positive gives a condition that can 
be written in the form 9i -t 9! -t 93 < 2(9192 -t 9293 + 939i)- By investigating what happens if the “<” 
is replaced by an equal sign and then solving for any one of the 9*1 this condition can then be translated 
into the one quoted in the main text above. 
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and T 2 on the 3-particle color singlet are not simply proportional to each other. [We are not 
sure how to argue mathematically that (A7) must also hold outside of the weak-coupling 
limit, but a violation would mean that there is some instability in the problem.] 

Eq. ( 5 . 34 ): (CI 4 , 0 ^ 2 , Cfg, 0) and (A^, At| 0) are the same except 

for how the states are implicitly normalized: e.g. (4.24) vs. (A_|_, A_| A'^, A'_) = — 

A'j^) — A'_). The different 5 functions of A± and of Cij are related by the Jacobian 

of the corresponding change of variables (5.28) and (5.33), giving the determinant factors 
in (5.34): (CI 4 , 0 ( 2 ! = | det ay|“^(A^, A^| and |C 44 ,C 23 ) = |A)[_, A)l)| det ay|“h Note that 
we get I detOyl”^ (and similarly | detay|“^) on the right-hand side because there is a factor 
of I detOyl^^A associated with each of the two dimensions of the transverse plane. 

Eq. (7.12): Consider the short-time {t —)■ 0) divergence of any matrix element 
(a;i|A(t) i?( 0 )|a; 2 ), where the \uji) are energy eigenstates with energies Ui. Insert a complete 
set of states to write the matrix element as X]E('^i|A(t)|i 7 )(i 7 |i?( 0 )|a; 2 )- The short-time 
divergence will come from intermediate states with arbitrarily high E, and in particular 
E 3> uJi,U 2 . Fourier transforming from t to a; gives u = E — wi, and so we see that the 
short-term divergence is associated with purely positive frequencies u. Now return to the 
original correlator (a;i|A(t) i?( 0 )|a; 2 ) in t space, and consider the Fourier transform integral 
f dt {ui\A(t) i?( 0 )|ci; 2 ) taking it to u space. For ca < 0, the contour may be closed in the 
lower-half complex t plane. In order for a short-time divergence of the integrand at t = 0 
such as t~'^ not to contribute when a; < 0 , the prescription for that divergence should then 
be {t — which is the prescription of (7.12). 

Eq. (9.9): Following ref. [ 6 ], for example, the last equality of (9.9) is obtained by doing 
the y integral first after recasting the range (9.3) as 


< y < Aty 


(A9) 


The result is 

r dy d{At) r d{At) f ^ \ 

J yAt - J ’ 


(AlO) 


with L dehned parametrically by (9.2) in our case. Then take the parametric range of the 
At integration to be io to L. 


Eq. (9.11): In both the cases of A{x, y) and A{l—x—y, y) in (8.1), there are some inter¬ 
esting cancellations between the pole and the integral pieces in (8.2). Individually, both of 
these pieces scale with y as 

0{y^/‘^{a+(3)) (All) 

for small y (and hxed x), up to logarithms. When the two are added together as in (8.2), 
there is a cancellation among the leading contribution from a and (3 to give 


= o(y/2Q,) + o(y'2/3) + o(y/s.)- (A12) 

dx dy 


There is no similar cancellation for the terms involving 7 . In the small y limit, {a, (3, 7 ) scale 
with y as 

a ^ (3 ^ y , 7 ~l/^. (A13) 

In consequence, both (All) and (A12) scale as y~^^‘^, and so the cancellation of the leading 
a+(3 term does not completely cancel the leading behavior. 
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Appendix B: Details on reduction of states from N to N—2 particles 


1. Normalization of projected states 


We first turn to the normalization factors, such as Af in the dehnition (4.21) of \{Bij}), 
necessary to dehne our projection from N particle states to effective N—2 particle states. 
One of our constraints is that YliPi — 0) discussion we will encounter corre¬ 

sponding p-space (5-functions evaluated at zero argument: 




EPi=0 



(Bl) 


where V± is the 2-volume (area) of the transverse spatial directions. Final results will not 
depend on V±, which is taken to be inhnite. Our other constraint is that = 0, and it 


will be convenient to give a name V± to the analogous h-space (5-function evaluated at zero 
argument: 






Y,Xibi=0 


(B2) 


Final results will not depend on V±, which is also inhnite. 

The projection V of an iV-particle state \b) onto the subspace with YhiPi — 0 
Y^i^ibi = 0 is 


V\bi, ■ ■ ■ ,b]\f) — 6j2xibi,oy±^ [ |hi-l-Ah, • • ■ , 6jv+^h). (B3) 

JAb 

(Here, Sj2xibo,o is a Kronecker 6 rather than Dirac (5-function and gives 1, rather than V±, 
when = 0.) One may verify that V‘^\b) = V\b). Because V\b) above depends only 

on the differences bi — bj of transverse positions, it is a state that can be characterized 
as depending only on the values of Bij = [bi — bj)/{xi+Xj). Except for normalization 
convention, it is the state \{BY) that we dehned in (4.21). 

The normalization of P|h) is 

(b), • • ■ , b'^\V\bi, ■■■ M = VY [ {b[, • • ■ , b'j^\b,+Ab, • • ■ , b^+Ab) 

JAb' 

N 

Ai i.i 

N 

= (B4) 

i=2 

where bij = bi — bj and we have assumed that bi and b' have been chosen with Yi ^ibi = 
0 = Yi^ib'i- These last constraints, together with Yi^i = 0) imply that 

N 

^a;j(b'i - bji) = 0, (B5) 

i=2 

and so the N—1 (5-functions left in the hnal expression in (B4) are not independent. (B4) is 
therefore equivalent to N—2 independent (5-functions times a factor of (5*^^i(0). In order to 
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relate this factor to the normalization factors we have dehned, we must take care how we 
write things. Ignore (B5) for a moment and rewrite the i=2 ^-function in (B4) as 


- b2l) = 


N 


N 


- bii) - bii] 


= Xo 


'-i=2 

N 


i=3 


N 


- bi) -^Xi{b[^ - ba] 


'i=l i=3 

In the presence of the other ^-functions, this can be replaced by 

- b 2 i) xj ( 5 ^^) ( ^ Xi{b[ - bi 

^i=l 

Given the constraint (B5), this is 


(B6) 


(B7) 


(B8) 


Combining with (B4) gives 

(6;,... , 6V|P|6i,... , 6„) = ifU. J]i(2)(6.i - b'j 


T-r ,, ,, X x'^Vi 


N 


2^± 

V, 


J](I.+I,)-V(S., -S'l). 

(B9) 


2=3 2=3 

In order to get rid of the nuisance factors of V± and V±, we hnd it useful to dehne 

'G, \V 2 




Vx 


) 


V\bi, • • ■ , hAf), 


(BIO) 


so that 


N N 

2=3 2=3 

Combining (B3) and (BIO) determines the normalization factor M of (4.21) in the main text 
to be _ 

u = (1414)"^/^ (B12) 

Eq. (Bll) with N=2 and iV=3 directly gives (4.22c) and (4.22b) for the normalizations 
in those cases. For iV=4, it gives 

= xl{x^ + x^)-\x, + x,)-^ <5(2)(S3i - s'^) 42)(B4 i - B^i). (B13) 

One may then use the transformations (5.14) to compute the Jacobian for trading variables 
(.B 31 , Bii) in the J-functions for (C 34 , C 12 ), remembering that we typically rename the Bij 
as Cij in the iV=4 case. The result of this change of variables is (4.22a). 
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2. Cancellation of normalization factors 


As discussed surrounding eqs. (4.26-4.28), the final state in the iV-particle language of 
El® El is 

lend) = f |hf, hf) = (B14) 

Jb( 

Now consider the initial state in this language. We have not yet picked a precise direction 
for the ; 2 -axis in our calculation—much of the formalism has been designed precisely to avoid 
that. For simplicity of presentation in the present argument, pick it for now to be exactly in 
the direction of the initial parton at the initial time. Then the system starts as |Pi=0)(pi=0|, 
which in the language of El (g) El is 

\staTt) =V];^\pi= 0 ,P 2 = 0 ) = f \bi,b 2 ). (B15) 

Jb\,b2 

Here, the factor of is the initial state normalization factor, just as described for (2.1) 
in appendix A. 

Because the projection onto the subspace of J2iPi — 0 Yhi^ibi = 0 is preserved by 
time evolution and splitting, and because the final state |end) lies in this subspace, no harm 
is done if we also project the initial state onto this subspace. That is, 

(end| ■ ■ • jstart) = (end|P- ■ ■ jstart) = (end| ■ ■ -Plstart), (B16) 

and so we may replace 


jstart) —> Pjstart) = Vj_ V / |hi, b 2 ). 

J 61,62 


(BIT) 


Using (B3) and (B2), this is 


jstart)—)■ Uj_ / 6 xibi+x 2 b 2 ,o\bi+^b,b 2 +Ab) 
J bi, 62 ,A 6 


= ur V,-^ 


(5(xibi + ^ 2 ^ 2 ) |hi+Ah, h2T^h). (B18) 


J 61,625^6 

For the initial particle, {xi,X 2 ) = (—1,1), and the above evaluates to 

Istart) ^ (U^fA)-' [ Kb,) = (U^Ua)"'jend) = =K\), 

Jbi 

where the last equality uses (B12). Combining (B14) and (B19), 

(end| ■ • ■ jstart) = (j ■ ■ ■ j). 

All the nuisance factors of Vj_ and Vj_ have canceled, as promised. 


(B19) 


(B20) 
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Z> 3,^3 -64,^4 

^^63,^3 ^^63,^3 

b2,X2 -^b2,X2 bixi -^b2,X2 

blx [ -61,61,^1’- 61,^1 

(a) (b) 

FIG. 28: The notation used in (a) appendix B3 and (b) appendix B4 to label different particle 
states immediately before and after a splitting, (a) is similar to fig. 13, but it was convenient to 
introduce slightly more general notation here. The use of primes {') here is completely unrelated 
to their use in (6.1). 


3. {B\6H\) matrix elements 

Here we work out the formula for the matrix element {B\SH\). For simplicity and deh- 
niteness, we will focus on the case where the initial particle in the amplitude splits into two 
daughters, using the labeling conventions of fig. 28a. The starting point is the SH matrix 
element in the amplitude, written more conventionally in terms of the individual particles 
in the Hilbert space H (rather than IHI0 H), given by (4.29). In p-space, this is 

{p2,P3\dH\p’2) = gT2'^23 ■ P 23 = gT2'^23 ' {X 3 P 2 - X 2 P 3 ), (B21) 

which in b-space becomes 

{b 2 , b,\5H\b',) = -igT 2'^23 ■ - 2:2 V53) [s^^\b 2 - b') 5^^\b2 - bg)]. (B22) 

Now turn to ]H[(8)]H[ notation and our effective N—2 particle description. Use the dehnition 
(BIO) and (B3) of the projected states to write 


{B\6H\) = ^{1,2,3\VSHV\1',2') 
V± 


^{1,2,3\SHV\1',2') 

V± 


= (bi,b2,b3|(5iT|b;+Ab,b'2+Ab) (B23) 

JAb 


(where we have suppressed writing the Kronecker ^-functions—the initial and hnal positions 
are to be understood as both satisfying the constraint Yhi^ibi = 0). Using the amplitude 
matrix element (B22), this is 


{B\5H\) = ST2'^23 ■ (a^sVf,, -a;2V53) [ 5^^\b,-b[-Ab) 5^‘^\b2-b',-Ab) 5^‘^\b2-b,). 
V± JAb 

(B24) 

At this point, the simplest way to proceed that will most easily generalize to other SH 
matrix elements is to note that (B24) can be rewritten, using (B3), in terms of the N=2 
state normalization as 

{B\5H\) = • {x3Vt,, - X2Vt,)[{bl,b2\r\b[,b',)^,,^,J^^\b2-b3)], (B25) 
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where the subscripts {x'i,x' 2 ) on the matrix element show the two Xi values that are to be 
used for both the bra and the ket.^^ Now we can use the N=2 version of the normalization 
formula (B9) to get 

{B\5H\) = -%g{x'2fr2'^2z ■ (xsVb, - 6^^\b2-h). (B26) 

The action of the combination 0 : 3 V 52 — X 2 ^on a function of h 2 —h^ is the same as the 
action of Vb, where B = S 23 = {b 2 —b^)/{x 2 +x^) is the same 3-particle B we dehned in 
(2.29). Also, we can absorb the overall factor of {x 2 )^={x 2 +X 3 )^ into 5^‘^\b2—b3) to get 
with result 

{B\SH\) = -igr2'^2z ■ Vb 5^^\B). (B27) 

This is eq. (4.13) of the main text. 


4. (C41, C 23 |< 5 iT|S) matrix elements 


The analysis proceeds similarly for 5H matrix elements between A^=3 and A^=4 particle 
states, where we label particles as in fig. 28b. The analog of (B25) is 


{{C,,}\5H\B') = 


_igV± 

Vx 


r 


2'^23 


■ {x 3 Vb,-X 2 Vb,) [(hi, b 2 , b4\V\b[, b'2, h'3).',4 S^^\b 2 -b 3 )_. 


(B28) 

Noting that ^4 is the third hj in the bra in (B28), the N=3 version of the normalization 
formula (B9) then gives 


{{C,,}\6H\B') = -tg{x'2)\x', -t- x[)-^ 6^^\C^,-B') • {xs'^b, - X2Vb,) 6^^\b2-b3). 

(B29) 

The rest follows as before. Noting that x^ + x[ = Xi + xi in the notation of fig. 28b, we can 
write the result as 


{{C,,]\5H\B') = -tg{x^ + x,)-H^^\C,3-B')r2'^23-^^^^\C23). (B30) 

Using an appropriate permutation of the dehnition (4.23), 

IC’ 41 , C 23 ) = 1x4-t-Xil (B31) 

we have 

(C41,C23|57/|S') = -^(/|x4 + Xi|-l<5(2)(C4l-B')r2^^23- V<5(2)(C23), (B32) 

which is equivalent to (4.15). The difference between this and the previous result (B27) for 
{B\5H\) is a factor of 

\xi + xi\-^5^'^\Cii-B') (B33) 

associated with the spectators. This is equivalent to the diagrammatic rule that we presented 
in the main text at the bottom of £g. 16. 


Readers may reasonably wonder why these Xi values instead of some other. The reason has to do with 
maintaining the consistency of our normalization convention (B2) for (5(^ Xibi) for both 3-particle and 2- 
particle states. In the above expression, 62 = ^3 and (from longitudinal momentum conservation) = xi 
and x '2 = X 2 + X 3 , which means that the 3-particle expression X\bi + X 2 b 2 + x^b^ for the final state is 
identical, in this context, to the 2-particle version x'^bi -|- x' 2 b 2 used for the final state in (B25). 
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Appendix C: Connecting our SH matrix elements to standard formulas 

In this appendix, we show how to connect the formulas given in section IV E to the usual 
textbook result 


{-p±, p±, zE\6H\0, E)rel 


2 


2 gVi 

2(1 - z) 


PU), 


(Cl) 


where the subscript “rel” indicates that the states are normalized using relativistic normal¬ 
ization. (This use of the notation ‘V” is unrelated to our 2 ; = 1—x—y elsewhere.) 

We first note that the Xi=l case of our splitting functions (4.34) correspond to the usual 
DGLAP splitting functions, with the even more usual helicity-averaged/summed result being 
given by 


(C2) 


For Xi different than 1, we can factor out the Xi dependence from (4.34) and rewrite these 
expressions in terms of the branching fractions 


Z = Zk = 


Xk 

\Xi\ 


and 


2 ; = = 


_ 


\Xi 


of the particular splitting. The relationship is 


P{Xi,Xj,Xk) = \Xi\P{-l,Zj,Zk) = \Xi\P{z) 


(C3) 


(C4) 


for each helicity case. 

Now consider the square of our matrix element (4.29), 


{Pj,Pk\SH\pi) 


2 


g^\Ti^jk ■ PjkY 


(C5) 


implicitly summed/averaged over color here. Using our dehnitions (4.30) and (4.31) of P 
and P, this is 


{Pj,Pk\5H\p,Y) 


2 


g^\e(±)-Pjk? 

Axjx'jxlE^ 


P{Xi,Xj,Xk) 


(C6) 


(where for brevity we again suppress helicity indices on P). For the choice of axis made in 
(Cl), where {pi,Pj,Pk) = (0, -p±,p±), we have Pjk = XkPj-XjPk = -{xk + Xj)p± = Xip±. 
Because Pjk is real, we have |e(±) ■ Pjkl"^ = Pfk- Then, converting (C6) to relativistic 
normalization by multiplying by {2Ei){2Ej){2Ek) = S\xiXjXk\E‘^, 


{Pj,Pk\5H\pi)^ei 


2 


‘^g‘^pi\xi 

\XjXk\ 


P{Xi,Xj,Xk). 


(C7) 


Using (C3) and (C4) then reproduces the standard result (Cl). 
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Appendix D: More details on divergences 


1. At—)-0 limit of xyyx 


Here, we will give some more details about how to take the small At limit of xyyx that 
gives rise to (5.46). For small At, cot(fiAt) and csc(fiAt) both become (At)“^ + 0(At). 
Then (5.42) becomes 


A, i; 

A, A, 
A, A, 


Y _ w _ 
^yy yy 

V - 7 - 
yy ^yy 



-IT 

y 


-1 


'' -IT -1 


(\M,\n, o\ 

V 0 V 

(\M,\n, o\ 

V 0 V 


At ^ 


-1T„-1 


«y 


+ 0(At). 


+ 0(At), 
+ 0(At), 


(Dla) 

(Dlb) 

(Die) 


Except for the MVL terms (coming from the evolution before and after the At interval), this 
is the result one would get in vacuum (ftj- = 0). Now use the simple formulas (5.16), (5.33), 
(5.35), which do not depend on details of the normal mode solutions, to get 


Ay A, 

A, A, 

A, A, 
A, A, 


Y - Y - 

_yy ^ yy 

Y - 7 - 

^ yy ^yy 


_iE{xi+Xi) fxiXi 0 \ f\Mi\Vt-, 0 \ 

At V -X 2 X 3 ) V 

iE{x3+X4) /X3X4 0 \ (\Mi\VLi 0 \ 

At V -X 1 X 2 ) \ 0 oy 

^ fX1X3X4 X4X2X4'\ 

At \X2X3X4 X1X2X3J ^ '' 


+ 0(At), 
+ 0(At), 


(D2a) 

(D2b) 

(D2c) 


Using the dehnitions (5.2a) and (5.4a) of Mi and Mf, the combination AyAy — AU that 
appears in the integrals (5.44) can be written in the form 


^-2 _ XiX2X3xlE^ 
^ ^ ^yy 


(f], sgn Mi + ftf sgn Mf) + O ((At)°). 


(D3) 


The fact that Ay and Ay vanish at the order of interest means that (5.45) simplihes to 


dV 

dxdy\ 


ClalM.Mi 

327rM2 


(-£l£2£3)r4) 


d(At) 

Jo TW 


X 


A yyYyyi^alo T j312 A I2^ A Ayy( q; a y a 27 )/f AO ((At)*^). 


(D4) 


Using (D2) and the integrals (5.44), expanding in At, and using = 0 to simplify, gives 
eq. (5.46) of the main text. 


2. Triple time coincidence associated with 1/At 

Here we justify more explicitly the claim in section VD that the hlf/At divergent terms 
in (5.46) for xyyx arise from the hrst three splitting times tx, ty, and ty simultaneously 
approaching each other. 
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In deriving results for xyyx, we integrated the first splitting time from —cx) up to ty 
in (5.8) to get 





Bx=0 


iMiBy 




(D5) 


In order to determine whether the region where is close to ty is important, let us instead 
integrate only from ty — At to ty for some At: 


rh 

J ty—Ar 


dtx Bx (-®y 5 I ^x) 


Bx =0 


iMiBy 


exp(|Mir2ii?y cot(r2iAr)). 


(D6a) 


Comparing (D5) with (D6a), we see that the only effect that the change of integration region 
has on our results is to replace 


hli —)■ r2®®(Ar) = —i sgn(Mi) fli cot(r2iAr) (D6b) 


in (5.42) and therefore in the divergent pieces of (5.46). This change has no effect at all on 
the Qf/At term in (5.46), and so the Qf/At term arises from t^ no more than At away from 
ty. This conclusion breaks down only if we make At so small that hi?® becomes so large 
that the small-At expansion made in (D2) breaks down. That happens when At < At. So 
the Qf/At term arises from the case where t^, ty, and ty all he within 0{At) of each other. 
A similar argument holds for the Qi/At divergence, with the replacement 

p+OO pty + Ar 

/ dt^^ dt^ (D7a) 

Jty Jty 

being equivalent to 

hlf —)■ r2f®(Ar) = —isgn(Mf) hlf cot(r2fAr). (D7b) 


3. ie prescription for double splitting (especially xyxy) 

In this section, we will more precisely justify the application of the ie prescription in 
section VIIB 2, especially as regards the xyxy calculation. As discussed above and in the 
main text, the 1/At divergences arise when three times become coincident, e.g. {ty,ty,ty). 
In order to sort things out, it will be helpful to know the explicit dependence on those three 
times as they approach each other. Unfortunately, we have already integrated over one of 
those times in our derivation of results like (5.46), (6.6), and (6.14). So let us step back and 
undo that integration. We could go back and rederive all our results (including doing all the 
B and Cij integrals) without integrating over the last (or hrst) time, but happily there is a 
trick that will allow us to extract the time dependence we want from the integrated results 
that we already have. 


a. xyyx 

We start by discussing xyyx, since this was the exemplar to which we related all other 
results. We will focus on the case where the last three times approach each other, giving rise 
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to the r2i/At divergence. In (D7), we saw how to integrate the last time over the interval 
(ty, ty + Ar) instead of the interval {ty, oo). If we take that result and then differentiate with 
respect to Ar, it is equivalent to never having integrated over tx in the first place. Using 
(D7b), this operation is equivalent to 


■ ■ ■ = —— / dt 

dAr 





(D8) 


with the understanding that Ar should be re-interpreted as tx—ty after all of these operations 
have been completed. We are interested in what happens when the last three times approach 
each other, which is the case where both Ar and At are small. So we may use the small Ar 
limit of (D7b), giving 


d 

dAr 



^ —2 sgn(Mf)/Ar 


(D9) 


Now carry out this procedure on (5.45) to recover the tx integrand. Then expand the 
integrand for simultaneously small At and small Ar, treating them as the same order. We 
hnd the result^® 


dr 


dV 

dx dy 

xyyx 

dx dy 


i^liSgn{Mi)xjx 2 xixi{xi+X 4 y{xs+X 4 )^ / d{At) / (i(Ar) 


Cial 


16n^ 

A{At+AT 


xyyx 

. 2 - 2 - 2 - 2 


X 


D3 


■ X2X4{xiXz-X2X4){oi+i3+'^) + 27 ) - (2a-h3/9-l-37) 


D 2 


D 2 


(DlOa) 


where At = ty — ty, At = tx — ty, and 

D = —X 2 X 4 At + {xi + £4) (£3 H- X 4 ) At. (DlOb) 

One may check that doing the Ar integral indeed recovers the Q-JAt divergence in (5.46). 
The details of (DlOa) will be mostly irrelevant to the present discussion; All we need to 
take away is that the only singularities of the integrand are associated with inverse powers 
of the denominator D dehned by (DlOb). D is proportional to the leading behavior of 
(AyAy — Ayy)(At)^Ar [with Of —)■ —isgn(Mf)/Ar]. And so one could have gotten what we 
will need for discussing ie prescriptions simply by looking at the structure of AyAy — AA 
rather than deriving (DlOa) in detail. 


3® Note that here [dT/dx dyY'^, which is the g —>■ 0 limit of dT/dx dy, does not correspond to all tl’s —>■ 0 
because has been replaced by —i sgn(Mf)/At, which does not vanish as g —>■ 0. 
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b. Permutations, and the Oi/At divergence of xyxy 

We may write (DIO) more explicitly in terms of the last three times by writing (DlOb) 


D = -{Xi+ £ 4 ) (£3 + Xi)ty + XiXsty - X2XAt^ 

= XlXsity - ty) - X2X4it5i - ty) (Dll) 

and rewriting the integration f d{At) d^Ar) as J dty dty dt^ divided by a factor of total time 
(and restricted to the appropriate ordering of the three times). In this form, the result turns 
out to describe the Qi/At divergences of not only xyyx but also of xyyx and xyxy as well. 
As a check, one may integrate the result over the last time in the xyyx and xyxy cases 
(tx and ty) appropriately and obtain ( 6 . 6 ) and (6.14) with At = ty — ty and At = ty — ty 
respectively.^® 

Now let’s consider ie prescriptions in the denominator D of (Dll). The times ty and t^ in 
the conjugate amplitude should have small negative imaginary parts compared to the time 
ty in the amplitude, as explained in section VII. For completeness, however, we should also 
consider the ie prescriptions of ty and t^ relative to each other. Two operators associated 
with times ti and ^2 both in the amplitude should be time-ordered, and so 

t 2 - ti ^ t 2 - ti - ie Sgn(t 2 -ti). (D12) 

Two operators with times ti and t 2 in the conjugate amplitude should be anti-time-ordered, 
corresponding to the conjugate of (D12), 


t 2 - ti ^ t 2 - ti + iesgn{t 2 -ti). (D13) 

A quick way to remember all of the ie prescriptions is to think of a Schwinger-Keldysh contour 
that runs from t = —00 to t = -|-cx) for times in the amplitude and then turns around and 
runs back again from t = -|-cxd to t = —00 for the times in the conjugate amplitude. The 
rule is that any time ^2 that appears further along this contour than another time ti should 
have a more negative (infinitesimal) imaginary part. 

So, for instance, we can implement all of the relative imaginary parts for the relevant 
times in xyxy (for which ty < t^ < ty) by 

(ty, tx, ty) y (ty f (^l P 62), tx ^*^2) ^y) • (^1^) 

Using the explicit values (4.7) of the £*, (Dll) can be put in the form 

D = (1 - X - y)(tx - ty) + (1 - a;)(l - y)(ty - U), (D15) 

for which the prescription (D14) gives 

D ^ (1 - X - y)(tx - ty + iei) -h (1 - a;)(l - y)(ty -tx + * 62 ) = D + ie. (D16) 


On a related note, the formula (Dll) for D is invariant, up to an overall sign, under (i) Xi —>• xj with 
(ty,ty,tx) —>■ (<y,ty,tx) sud uuder (ii) Xi —>■ ici with {ty,ty,tx) —>■ {ty,tx,ty). The / dty dty dtx integral 
representation of (DIO) is invariant under the same transformations with additionally (i) {a,/3,"f) —>■ 
(d, a, 7 ) and (ii) {a, (3, 7 ) —>■ ( 7 , a, (3) respectively, as in section VI. 
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The important point here is that we could get the same prescription by instead writing 

D ^ {I - X - y){t^-ty + ie) + (1 - a;)(l - y){ty - tx), (D17) 

which is At —)■ At+ie for the At = t^—ty relevant to xyxy. This agrees with the prescription 
guessed in section VIIB 2 of the main text for the Qi/At divergences of xyxy. [In this case, 
it happens to also be the same as the prescription one would get from simply applying (D13) 
in isolation to At = tx — ty. However, that is an accident, as we will see in a moment.] 

One may similarly check the ie prescriptions used for all of the other Oj/At terms in 
section VIIB 2. 


c. Of/At divergence of xyxy 

Now consider the Of/At divergence of xyxy, which arises from the first three times 


(tx,ty,tx) approaching each other. The corresponding value of 77 oc {XyXy 
[now with Oj —)■ —isgn(Mi)/Ar] is 

-A2-)(At)2Ar 

D = —X 2 X 4 At + (yXi + 74 ) (xs + X 4 ) At, 

where At = tx — ty and At = ty — tx. This can be written explicitly as 

(D18) 

D = {1- a;)(l - y){ty - Q + {1 - x - |/)(tx - ty). 

The relative ie prescriptions for tx < ty < tx in xyxy can be organized as 

(D19) 

(tx, ty, tx) y (tx, ty i(ei T *^ 2 )) tx i^ 2 ') ■ 

It is convenient to reorganize (D19) as 

(D20) 

D = (1 - x){l - y){tx - tx) + xy{ty - tx), 

from which is it particularly easy to read off that (D20) gives 

(D21) 

D ^ D — ie. 

(D22) 


This is the same prescription as just taking At —)■ At —te in (D19), which matches the guess 
that we made back in section VIIB 2 of the main text. [Note that this is not the prescription 
we would have gotten by applying (D13) in isolation to At = tx — ty.] 

One may similarly check all of the other 1/At divergences of section VIIB 2. 


Appendix E: Relating different crossed interference contributions 

In this appendix, we show more explicitly how to relate the xyyx and xyxy contributions 
to the xyxy contribution. 
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1. xyyx 


a. The basic discussion 


For xyyx, the analog of the xyyx formula (4.16) is 


dl 

dx dy 


xyyx 




E 




B ^=0 


X ^ l-x-y,x) ■ 

X ^ ^-y,y) ■ ^cl){Bh^h,,hy{l-x l-x-y,y) ■ V 

X Phi^ft,hx(l -t l-X,x) • V B^{B^ ,ty\B'",ty,) 

B^=0 


(El) 


[where we have used the notation of (4.37) rather than (4.36) for the splitting functions B] 
Substituting (4.38) for (4.37) and then using (5.9), the analog of (5.10) is 


dV 



dx dy 

xyyx 

Stt'^E'^ \x'i + x'^Wx', -I- 0:4 


Bl Bl 


exp(-i|Mf|fif(i?")' - \\M,\n,{B^f) 

By,By 

X Cl,, ^t\Cl, Cl„ 0 ) 


"" Jo JB^,By {Byy {Bff 


C'i2=o=C'?3; cl.=By-, cl.=By 


(E2) 


with At = ty — ty. The Cij and 4-particle propagator are now implicitly dehned in terms of 
Xi = x\ instead of Xi = The derivation of (5.39) is symmetric in such a way that 


(CI„ Cj's, Atici, C\„ 0) = (Ci, c;„ Af lCl,C^, 0). 


(E3) 


Using this relation, the xyyx result (E2) is superhcially the same as making the substitutions 
summarized in section VIA on the xyyx expression (5.10). We say superhcially only because 
there is a sign issue that must be addressed when replacing Xi —)■ x[ in the formula for 
(CI 4 , CI 2 , C 23 , 0 ), which we discuss next. 


b. A sign issue 

When the 4-particle Xi are the (4.7) appropriate for xyyx, then there are a number of 
important and related properties of the 4-particle normal mode formulas in section VB2: 

• = + + + \/A is always positive, so that the normal mode frequencies 

12+ and hi- given by (5.21) are both proportional to 

• XiX 3 / 2 N^E and XiX 3 / 2 N_E are both positive, so that the normal modes (C* 34 ,C*i 2 ) 
and (C^, 0 ^ 2 ) of (5.24) are normalized to be real. 
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• As a result, the normal mode degrees of freedom A± defined by (5.26), which were 
implicitly assumed to be real, are indeed real. 

In contrast, when the 4-particle Xi are the x' (6.1) appropriate for xyyx, the situation 
changes: 

• A+ dehned above is still positive but A_ is negative, so that 12+ oc but hi. oc \/+*. 

• X 1 X 3 / 2 N+E is positive but xiXs/ 2 N_E is negative, so that the normal mode (C^, Cj^) 
has been normalized to be real by (5.24), but the normal mode (^^,^ 12 ) has been 
normalized to be imaginary. 

• As a result, the dehnition (5.26) makes A_ imaginary instead of real. 

Imaginary values for (C'^,C'i 2 ) confuse the analysis that follows in section VB3, especially 
the normalization of the relation (5.34). The easiest way to clear this up is to keep to real 
normalizations for the normal modes: 



and 

(A+,A„)“-^(A+,-^A.). 


This means that the Lagrangian (5.29) for A± is then normalized as 


(E4) 

(E5) 


± 


(^new)2 


nl 






(E6) 


Now the A+®® are both real, but the Lagrangian for AA®" has an unusual overall sign. To 
understand the effect of the sign, it is useful to look at the corresponding Hamiltonian 


« = E ± [5 + 5 

± 


(E7) 


In time evolution exp{—iHt), negating a Hamiltonian has the same effect as negating time 
t. So the A_ part of the time evolution (5.30) changes from 

(A_,t|A'_,0) =- —— exp^i[|(A^-|-A'^)r2_ cot(r2_t) —csc(r2_f)] j (E8) 

to its conjugate 
(A_,t|A'_,0)“" = 

- exp(^-i[i(A2 + cot(H_t) - (A_ ■ A'_)“"'H_ csc(H_t)]). (E9) 

The one change here that matters is the appearance of an overall minus in the prefactor, 
which we will address in a moment. The other change is to the sign of the argument of 
the exponential, which we now show is inconsequential. Combined with (E4), its effect is to 
change expressions of the form 

cot(H At) and csc(H At) (ElO) 
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in the propagator (5.39) to 

flneJ (^ cot(ll At) and ^ ) Qcsc{QAt) (Ell) 

where, based on (5.28), (5.33) and (E4), 

anew = a(^\)- (E12) 

Using the fact that is a diagonal matrix (5.40), the new expressions (Ell) are the same 
as the old (ElO). 

The one important change is the overall minus sign in the prefactor of (E9), but we have 
already accounted for it in the main text when we decided in the formula (5.39) for the iV=4 
propagator to rewrite the overall factor \xiX 2 X 2 ,X 4 \ as —X 1 X 2 XSX 4 . For xyyx, for which the 
Xi are Xi (4.7), this had no effect because 

—XiX 2 XzXi = \xiX 2 XzXi\. (E13) 

For xyyx, however, the Xi are x[ (6.1), and 

—x\x' 2 x'^x'^ = — \x\x' 2 x'^x'^\-, (E14) 

so using —xiX 2 X^Xi instead of \xiX 2 X^X 4 \ introduces precisely the additional overall sign that 
we need to account for in the prefactor of (E9). 


2. xyxy 


A similar analysis applies to xyxy. The analog of the xyyx formula (4.16) is 
dl 

dx dy 


J xyxy 


= d-^^ tr(T“TKE^)^|5^i + 54|-1 x3 + 54 


-1 


E 

h,,,hy,h^,h,h 

X Vh^h,,hy{l-x l-x-y,y) ■ VBy{B^,tj\B'',U) 

By=o 

X {Vh,^h,hy{^ ^ 1-1/,I/) ■ '^c!j{'PUh.,hA^-y ^ l-x-y,x) ■ VCI3) 

{cl,cl,Mcl,c!2,h) 


X Vhi^h,h,,{l l-x,x) ■ VBx(B^,ty|B^, 4 ) 
and the analog of (5.10) is 




(El 5) 


dr 

\_dx dy 




xyxy 


X 




d{At) 


\ ‘^4||‘^3 ' ^41 

" exp(-||Mf|fif(i?^)'-i|Mi|fii(i?")2) 


iQ JBy,B^ 

^ Qa, At|Cf4, CI 2 , 0) 

12 


(E16) 

ICf2=0=Cf3; Cl^=By-, Cfi=B^ 

with At = — This result is related to the other diagrams by the transformations given 

in section VIB. 

There is also exactly the same sign issue as the one discussed above for xyyx, which is 
accounted for in the same way by our choice of using —X 1 X 2 XSX 4 instead of \xiX 2 X 2 ,X 4 \ in 
(5.39). 
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Appendix F: More details on comparison with Refs. [6, 8] 

Here we give a few details on making contact with Blaizot and Mehtar-Tani [6] and Wu 
[8], which explored energy loss in the limiting case y ^ x 1. We will focus mostly on the 
comparison with Wu, which uses language closest to our own. 


1. The double log region (9.3) 

In Blaizot and Mehtar-Tani [6] (50),^° the double log arises, in their notation, from 
qr Using their dehnition^^ r = ca'/g^, these conditions may be rewritten as 


u 

^ < r < 

pZ 



(FI) 


The relevant is identihed as ~ y/uJq after their eq. (77). Their ca, ca', and r correspond to 
our xE, yE, and At. With this change of notation, (FI) above becomes our (9.3). 

In Wu [8], the second inequality of our (9.3) is given explicitly in his (38). The hrst 
inequality is used implicitly in the expansion of his (43) into (44), which assumes [in addition 
to the explicit conditions of Wu’s (38)] that uj 2 x\/t 2 1 (in Wu’s notation). Given that 
x\ ~ ^/qL, the two conditions can be written as 


U2 

^ < t2 < 
qL 



(F2) 


In our notation, 1 x 2 and t 2 are yE and At. In our thick-media approximation, the relevant L 
is the formation length ~ ^yxE/q of the x gluon, which then makes (F2) above equivalent 
to our (9.3). 

In our own work, the second inequality of (9.3) corresponds to max(r2_|_ At, At) 1 
and so corresponds to the small-At expansions of (D2) [but one must include more terms 
of the expansion than shown there to obtain (9.4)]. In contrast, the hrst inequality of (9.3) 
guarantees that XyXy — XU is not given by its small-At expansion (D3) in the cases of 

{xi,X 2 , X 3 , X 4 ) = (-1, y, 1-z-y, z) = (-1, y, x, 1-x-y) (F3) 

for zyyz [given hj x ^ z = 1—x—y in (4.7)] and 

(xi,X 2 , xs, X 4 ) = {-(l-y), -y, 1-z, z) = (-(1-y), -y, x+y, 1-x-y) (F4) 

for zyyz [similarly from (6.1)]. The different behavior of XyXy—Xyy is due to an approximate 
cancellation of the leading terms when y ^ x ^ 1. 


Similar to our convention for Wu (our footnote 31), we number equations from Blaizot and Mehtar-Tani 
[6] according to the arXiv version of their paper (version 2). 

They write t = uj/q^ just before their (50), but the relevant frequency in their later discussion of radiation 
in section 4 is what they call w' there. 
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Wu 

y < X <C 1 limit of our zyyz 

z 

t 

UJi 

xE 

102 

yE 

{Ki,K2) 

(U_,U+) 

(Hi, U 2 ) 

(^ U_, 0+) when y <C x <C 1 

ix,B) 

(^34,^24,) = ((a: 3 -|-X 4 )C 34 , —X 1 C 12 -I-X 3 C 34 ) ~ ( 0 * 34 , 012 ) 

my^x 



A+ 

G(3) 

{B,t\B’,0) 


(A+,A_,t|AV,A'_,0) 


TABLE I: Some translations between Wu [8] and the y < x <C 1 limit of the zyyz case of this 
paper. 


2. y < X 1 formulas of Ref. [8] 


Let’s focus on our zyyz, which corresponds to the second diagram in Wu’s eq. (30) for 
his Cl contribution to energy loss. Wu’s formula for our dV/dxdy, extracted from his (33), 
would be 


dl 

dx dy 


zyyz 


2uj\uJ2 


/ / d^X 2 d^x^ 

' Z\<Z2<ZZ<Z4, J 

£C 4 lG^^^(a;4, 2:4; ccg, Z3;ui) G^^\x2, Z2; Xi, Zi;ui) 

I B2=X2,B3=0 


1xi=0=X4 


X 


Vb2-^BsG^^Hxs,B,, Z:,-,X2,B2,Z2;u}i,u}2) 


(F5) 


where the notation on the left-hand side is ours and the notation on the right-hand side is his, 
except that the hrst integral sign is just our way of writing the J dz^ dz^ dz 2 dzi integration. 

and are given by Wu (14-16) and (20). A translation table between Wu’s notation 
and the y ^ x limit of our notation is given in table I. When comparing to our formulas 
in the main text, keep in mind that the 4-particle Xi for zyyz are given by (F3) and not 
(4.7). Using the translations of table I, one may verify that Wu’s formula (21) for relating 
{x, B) to {x, B) is equivalent to the y ^ x <^1 limit of our formula (5.27) [with (5.24) and 
(5.28)] for relating A± to (C 34 , C 12 ). 

With these translations, we hnd that the y ^ x limit for the zyyz version of (4.40), 
combined with (5.38), reproduces (F5). However, for y ~ a; 1, we hnd some differences. 
One difference is that Wu (20) is missing a factor of the Jacobean of the transformation 
between the variables {x,B) and {x,B) [as can be seen by taking the Z 3 Z 2 limit, in 
which all the Green functions become J-functions]. The corrected relation can be written in 
the form 


G(^H^3,B3 ,Z3;X2,B2,Z2;CJi,CJ2) - 

- G(* 3 , X2, Z3-Z2, mi, Ki) G{B 3 , B 2 , Z3-Z2,m2, K 2 ). (F6) 

a;ia;2 
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Another difference appears to be in the treatment of the triple gluon vertex in Wu (30). We 
hnd multiple terms, associated with the different helicity amplitudes in our section IV E and 
their combination into (4.38). The derivative associated with that vertex (our V^^y^) is also 
more complicated than Wu’s Finally, the two 3-particle evolution frequencies in (F5) 

should not both be Ui when x ^ y. Altogether, our result for y < a; 1 replaces (F5) by 


dl 

dx dy 


a^N.CAE^ 


J zyyz 


^nn^mm ^nm^nm ^nm^nm 
+ - 


2ufu2{uJi+U2y V UJ 1 +UJ 2 ^2 / Jzi<Z2<Z3<Zi 

1 xi=0=X4 

^4; X3, Z3; CUi) G^^\x 2, Z2] Xi, Zi, U1+U2) 


d‘^X 2 


V™ V” 

^ Xi ^ X4 


X {uiV^s^-U2VljV%G^^\xs,Bs,zs;x2,B2,Z2;u,,U2) 


B2=X2,Bs=0 


. (F7) 


Note that, in momentum space, the combination a;iVB 2 —‘-^ 2 '^X 2 above is (proportional to) 
simply one of our rotationally-invariant combinations Pij = XjPi — XiPj of the transverse 
momenta associated with that vertex. 


3. Problem with comparing the double log contribution of xyyx xyyx 

There is a subtle problem comparing (i) our results for the xyyx + xyyx contribution to 
the double log to (ii) corresponding pieces of refs. [6, 8]. (These contributions correspond 
respectively to the second diagram in Wu (31) [8] and the third diagram in Wu (32), and 
they correspond to and G 2 of Blaizot and Mehtar-Tani [6].) A simple indication of the 
problem can be found by comparing (i) our potential (5.19) applied to the x,y 1 limit 
of xyyx or xyyx to (ii) the potential of Wu (17). Consider xyyx, where our notational 
convention (4.7) was that ^2 is the position of the y gluon and h 4 is the position of the x 
gluon. In the x,y limit, the two harder particles stay extremely close together (relative 
to other transverse distance scales), corresponding to hi —)■ h^. Adopting Wu’s convention 
of using X for the separation of the x parton from the hard partons (our h 4 i ~ h 43 in this 
case), and B for the y parton from the hard partons (our 621 — 1)23 here), the corresponding 
limit of (5.19) is 

V ~-^(2a;2 + 2B2). (F8) 

8 

This does not agree with Wu’s 

V + + (F9) 

8 

The reason for the discrepancy is that Wu assumes that the two hardest particles form 
a color dipole with total color in the adjoint representation A. When those particles are 
gluons, however, they could in principle form a dipole in any color state contained in A(8)A, 
which contains many more possibilities than A itself (especially in the large W limit). This 
situation is depicted in Eg. 29a (which is drawn by distorting Eg. 11 into a style similar to 
Wu’s diagrams to depict relative transverse separations). In this picture, the initial singlet 
dipole Erst emits a gluon and must become an adjoint dipole to balance color. The emission 
of the second gluon in principle allows the dipole to become anything in A®A. However, all 
of these possibilities other than A itself will be suppressed by the small dipole size once we 
also add in virtual processes where the y boson connects to the other particle in the dipole. 
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xyyx 



(a) (b) 

FIG. 29: (a) A depiction of possible total dipole color at different times for xyyx for the case 
y < X <C 1. (b) A virtual correction to single bremsstrahlung that conspires to suppress the range 
of dipole color possibilities in energy loss calculations. 


zyyz 


-k. 


FIG. 30: Similar to fig. 29a but for zyyz (hg. 9). 


as in fig. 29b. Since in this paper we are jnst compnting dT/dx dy (for fixed x and y) and 
not energy loss, we have not inclnded virtnal processes, and so the other color possibilities 
are not snppressed. The color dynamics of onr calcnlation is therefore different than Wn’s, 
which makes it difhcnlt to directly compare partial resnlts. (The same holds for comparison 
to Blaizot and Mehtar-Tani.) 

For zyyz, however, all is well, as depicted in hg. 30: In the x,y 1 limit, the calcnlation 
of [dT/dx dy]zyyz is also restricted to an adjoint color dipole at intermediate times. On a 
related note, following the nnmbering convention of (F3), ^4 —)■ hi for <C 1, and onr x 
and y gluons correspond to 63 and 62 respectively. Wu’s x is then our 634 ~ 631 , and his B 
is our ^21 — ^ 24 . In this limit, our potential (5.19) then indeed agrees with Wu’s potential 
(F9). 
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